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Diffusion dans les couches limites stratifiées 


(continuação) 


Nous allons par la suite donner un certain 
nombre de résultats originaux montrant assez bien com- 
ment une action thermique destabilisanta est susceptible 
d'accélérer la processus de diffusion. Avant d'analyser 
ces résultats il convient de donner três brévement une 
idée des difficultés théoriques rencontrées dans I'étude 
d'un mécanisme qui se trouve situé en quelque sorte au 
coeur même du problême turbulent. Deux voles d'appro- 
che s'offreni à nous: soit une théorie macroscopique 
qui utilise toujours en fin de compte una analogie molé- 
culaire soit une théoria statistique qui cherche à pénétrer 
plus à fond le phénomêne mais se heurte d'embléa à la 
difficulté centrale du problême turbulent: son caractâre 
non linéaire. 


H.l L'analogie moléculaire 


Les premiers efforis pour constituer un modéle ana- 
Iytique de diffusion turbulente ont procédé par analogie 
avec le phénomêne déjã connu de diffusion moléculaire 
tel qu'il est exprimé par exemple dans une loi du type 
«conduction de chaleur» alors a — 1/:e, 


dc à [ de” 
= a 
dt 0x, 


z est le coefficient qui rêgle la ciffusion de C. Cette loi 
sous entend une loi du type Fourier oú le flux (à travers 
tout élêment de surface) de la quantité est liê à la pro- 
jection du gradient de C suivant la normale à cet élément, 
L'isotropie du miligu est également postulée ce qui 
éêvito de donner à a un caractêre tensoriel, 


La linéaritê du problême est conservée si a dépend 


des x, mais non de C. Si a est constant il est possible 
de mettre en évidence l'opérateur diffusion et d'écrire 
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d. MATHIEU 


ot OX 

Dans la cas des couches limites disons que «grosso 
modo» un processus de diffusion transversal est associé 
à un processus de convection longitudial. L'équation qu' 
contrôle |"évolution d'un «scalaire passif Cx» est: 


d o Fa per A ad o a 
C+u. = A | o (x) a Cc | 
] 


í a : ' | C 
+ U e (x:) du + eu lx: ci 
E * e EM O e 
| ] L 

avec UC=-—e, dc/0x, 


Le symbole de Kronecker ER fait ici apparaitre le 
caractôe supposé isotrope de la diffuslon molécu- 


laire tandis que le tenseu; E (x) souligne le carac- 
têre complexe du phénomêne turbulent; un flux dans la 
direction «jr» n'étant pas nécessalrement associá à une 
variation de C dans cette même direction. Cependant 
la processus qui associs au flux suivant j un gradient 
dans la même direction risque bien d'être le plus actif. Si 
Ion postule alors l'isotropie du mouvement turbulent 
lui-même 


et Ion a à nouveau: 


—T+U, — = - | a (x) + e (x) 


| 
Ox, 


mais l'on peut moins sommairement se contenter de ne 
retenir que les éléments diagonaux du tenseur Ei ils 
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sont souvent appelés coefficients d échange s;, es, es bien 
que |'expression coefficients principaux de diffusion 
turbulente paraitrait plus parlante. Le plus souvent on 
néglige aussi dans cette formulation la diffusion molé- 
culaire. 


Dans les cas qui nous intéressent on peut faire 
les hypothêses suivantes: 


- ja champ de concentration est stationnaire 
9C/t=0 
— [Vécoulement moyen est parallêle à l'axe 1 
U=U,=0 
- la diffusion longitudinale (c'est-ã-dire suivant l'axe 


1) est négligée vis-à-vis du transport convectif lié à la 
vitesse moyenne. L'équation prend alors la forme 


= BE i ac ê alo 
Elk === Es TT E <A Em. 
Ox, Axo | Oxta | OXa | Ox | 


Dans le cas d'un écoulement de couche limite sur 
paroi plane dans lequel une source ponctuelle de con- 
taminant situé au niveau de la paroi injecte le débit O 
les conditions aux limites sont les suivantes 


ie 


si x —> 0 C —+ o 
si = Cc 1 

Sc 
Pour x >0 et x 0 Es —> O 


Enfin y x, on doit avoir 


+: MI” 
— O | | U, lx) C (x) dx; dx; = Q 


. 


si U,, e: 6t e; sont des constantes, le systâme admet la 
solution suivant pour x) >> x” + xy 


0 | U, Xu Ra. | 
C (x) = / = GRE: 1 sos E - 
| Â x: tor Es | 


E sr Ka V Ea Ega 


La solution proposée s'est révélée en net désaccord 
avec |'expérience les mesures indiquant que la concen- 
tration maximale décroit en fait beaucoup plus rapide- 
ment que ne le prévoit le calcul. Certes, il est possible 
de donner à l'équation de diffusion un tour moins 
schématique surtout si l'on tient compte des movyens 
actuels mis à disposition pour le traitement numérique. 
Cependant les nombreuses tentatives qui ont eu lieu 
pour rendre plus réalistes les hypothêses de base n'ont 
pas donné les résultats escomptés. || n'est pas douteux 
que le cadre dans lequel on s'efforce de circonscrire 
la problême est quelque peu simpliste. L'introduction 
de la notion de longueur mélange n'améliore pas la situa- 
tion si même ella ne la détériore pas, En effet, si du 
coefficient de diffusion turbulente : on ne retient que 
sa dimension longeur vitesse on est tenté de le lier à une 


238 


dimension «l» et à une vitesse «u» toutes deux caracté- 
ristiques du champ fluctuant. On pose alors: 
e-IVu 

(u peut être la fluctuation de vitesse et | une longueur 
de corrélation). Mais avec la notion de longueur de 
mélange dont l'idés de base est empruntée à la théorie 
cinétique des gaz on interprête «l» comme la distance le 
long de laquelle la quantité transférable (quantité de 
mouvement, matiêre...) circule sans altération. Plaçons 
nous dans la cas d'une couche limite plane définie dans 
la plan 1. 3. Par exemple surgie du point d'altitude x; 
ladite quantité se réincorpore brusquement au milieu 
ambiant au point d'altitude x; + |. Ceci invite à lier 
la fluctuation u, au champ moyen U, ce qui se traduit 
en premiêre approximation par la relation: 


du, 


d xs 


Prandtl admet comme vraisemblable l|'hypothése 


a 


selon laquelle us — u; Si es — ls V u!, 
dU, 
EX 


il viant e, — Fa 


et la coefficiente de diffusion turbulente se trouve ainsi 
liá au gradient du champ moyen. Or, il est bien connu 
que tout champ turbulent diffuse même s'il est associé 
à un écoulement uniforme. La contradiction est évidente 
du reste, toute relation qui tente d'établir ponctuel- 
lement une relation entre les propriétés du champ turbu- 
leni et celles du champ moyen recêle quelques Ines de 
ces difficultés. 


Réflexions sur les probêmes de diffusion 


Avant d'aborder le point de vue statistique, il con- 
viegnt d'analyser sommairement les faiblesses de |'ana- 
logia moléculaire et de prendre conscience de la diffi- 
culté centrale du problême posé: 


|| est toujours sous entendu, lorsque Ion aborde un 
problême de physique mathématique décrit par une équa- 
tion cifférentielle, que les évolutions spatiales ou tem- 
porelles de la quantité traitée mettent en jeu des êéche- 
les de longueur et de temps infiniment plus grandes 
que celles nécessaires à la définition des coefficients 
figurant dans |'équation. Dans un problême de conduction 
de chaleur les caractéristiques thermiques du matériau 
traduisent des propriétés structurales définies à échelle 
microscopique, les temps de relaxation sont également 
três faibles; quant à l'équation elle ne cherche quã 
traduire une certaine évolution macroscopique de la 
température et son traitement ignore geênéralement 
tou: des aspects propres à la physique de la matiêre. 
Si l'information thermique est transmise instantanément 
en tout point, léloignement à la source perturbatrice 
ne faisant que régler l'intensité de cette information, 
cela suppose un mécanisme de transmission sousjacent 
dont la constante de temps est négligeable vis-à-vis 
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de celle du phénomêne principal. Cette situation privi- 
légiée me se retrouve absolument pas lorsque la diffusion 
est assurés par voie turbulente. Dans une couche limite 
cinématique les plus grosses structures turbulentes ont 
des dimensions qui sont à peine plus petites que celles 
de la couche elle même et leurs temps caractéristiques 
sont tout à fait comparablos à ceux qui interviennent 
dans les évolutions moyennes du systême*. On doit donc 
penser que toute tentative cherchant à prendre en compte 
ces effets au moyen d'un simple coefficient « est vouée 
à un échec. Le fait de reconnaitre à r un caractêre ten- 
soriel suppose une anisotropie fonciêre du milieu mais 
n'améliora pas la situation du point de vua qui nous 
occupe. Si les tentatives de liaison entre champ turbu- 
lent et champ movyen se sont avérées en rêgle générale 
malheureuses cela vient du fait qu'elles recélent d'un 
certain point de vue les mêmes insuffisances. Les bouf- 
fées turbulentes ne s'adaptent pas instantanément au 
champ moyen local bien qu'elles le modâlent à travers 


les corrélations u; U,; car dans ce chaos une certaine 
spécifité demeure toutes les strutures ne contribuant pas 
également à une caractéristique turbulente donnée. || est 
maintenan bien connu que la contrainte tangentielle 
Tig Nest pas nécessairement nulle en tout point oú 
à U,/ôx l'est. J. Mathieu [20] a montré qu'une telle 
situation apparaissait dans le cas d'un jet pariétal infir- 
mant tout caractêre général à des relations du types 


3U, 


O Xs 


Une explication recevable a été donnée: [21] les 
structures turbulentes nées dans la parties extérieurs du 
jet inscrivent en partie leur histoire dans la région parié- 


* Comme nous l'avons déjãá maints fois signalé, il est de ce 


tale, tandis que s opêre à la fois leur restructuration 
contrôlés par l'effet cumulatif des déformations, leur 
mutation par émisttement et disparition par frottement 
visquex. L'évolution du coefficiant da corrélation du fait 
du gradient de vitesse a été extraite des travaux expéri- 
mentaux de Champagne et Corrsin [3] et de I'équation 
de Craya [7] traités par Courseau et Loiseau [6] dars 
la présentation donnée par J. Mathieu [21]. 


Dans le cas du jet pariétal, I'échec d'une relation 
simple du type Boussinesq est entre autre imputable au 
phénomêne de convection mais on connait des cas tel 
que celui relaté par Hanjalic et Launder [13] ou les 
termes dits de diffusion turbulente sont les seuls res- 
ponsables. Ceu auteurs opérent en effet dans un canal 
plan à parois fortement assymétriques** et constatent 
également un décalage entre point de vitesse maximum 
et point de tension nulle. Du fait que dans ce cas les 
lignes de courant moyennes sont parallêles à l'axe du 
tuyau, le régime étant établi on ne peut qu'évoquer 
influence des termes dits de diffusion turbulente pour 
saisir le mécanisme en cause, du reste se sont eux qui 
seront appelés à jouer le rôle essentiel dans le processus 
de diffusion qui nous occupe. 


La convection agit évidemment sur toute quantité 
turbulente mais la décomposition commode, bien qu'arbi- 
traire, du champ global en champ moyen et en champ 
turbulent*** nous amêne à distinguer la convection lige 
à lécoulement moyen (la seule que nous ayons prise 
en charge dans le jet pariétal) de cella lie à I'écoule- 
ment (cas du canal assymétrique). C'est ainsi que nous 
aurions dú écrire l'équation relativa à l|'énergis ciné- 
tique turbulente 


à a à é 0 ou, 
8 O Mid Mega MA Bird 
1 Ó | , u; = Ea 
= e dx; [75 | — “u; FA ' E 9 u; 65; 


le lar membre mettant en relisf non plus l'évolution lors- 
que Ion suit I'écoulement moyen mais bien cette même 
évolution pour une structure véhiculée à la fois par le 
champ convectif moyen et par le champ turbulent. L'exis- 
tenca d'un terme convectif turbulem lié au caractêre non 
linéaire de I'équation de départ rend compte en défini- 
tive du parcours aléatoire de la particule («random 


fait impossible que l|'ensembla des structures turbulentes se 


trouvent en état d'áquilibre ([«quasi équilibrium») dans un. systême de variables réduites constrult à partir d'échelles de longueur et 
temps définies à partir de cet écoulement moyen. Certaines structures turbulentes contrôlent en effet Ievolution de ce dernier et 
leur temps caractéristique est tout à fait comparable à celui relevé sur lo champ moyen. C'est ainsi que Sunyach et Mathieu [33] 
trouvent un profil intermittent qui transcrit dans les variables réduites classiques ne cesse d'évoluer traduisant une sorte «d'émiettement 
croissant du modêle. Cette évolution est du reste plus marquée dans le cas oú il s'agit d'écoulements movyens fortement évolutHs 
(jets). Même dans le cas d'une turbulence engendrée par grille G. Comte-Bellor et S. Corrsin [4] déácêlent des structures turbu- 
lentes présentant un caractêre pseudo-permanent dont il convient de tenir compte si Ion veut respecter les lois de décroissance 


observões. 


** Du falt de la non linéarité les deux modéles ne cessent de s'interactionner, c'est ainsi que les corrélations u: 


VM s“intro- 


duisent dans |'historis de I'écoulement moyen, celles d'ordre 3 apparaissant dês que Ion cherche à fixer I'évolution des précódentes. 


*+* Una des parois est lisse et I'autre fortement rugueuse, 
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paths). C'est dire en quelque sorte que la non linéarité 
se manifeste entre autre dans ces équations traduisant 
Hávolution des propriétés movennes du champ fluctuant 
à travers ce cheminement cahotique qui est |'essence 
même de tout problême de diffusion. Bien entendu il 
faut s'attendre à ce que la paramêtre pression soit égale- 
ment affecté par cette non linéarité puisque champ de 
pression et champ de vitesse sont intimenent liés. Dans 
la cas d'un écoulement à gradient de vitesse constant 
Craya [7] montre comment est étable cette liaison 
donnant dans l'espace spectrale une relation ou seules 
figurent des corrélations de vitesse (double et triple). 
A priori cependant aucune raison physique ne milite 
en faveur d'ure connection aussi étroite que cells pos- 
tulés dans la plupart des méthodes de calcul de la 
couche limite ou l'on écrit par exemple: (la terme dit 


de diffusion turbulente est associé au gradient q). 


= à 7 
Sm: Rs 
Us mes — ==g . 
Et Y À xs 


seul l'embarras causé par la présence de tels termes 
excuse une formulation aussi discutable. La distribution 


spatiale par exemple de I'énergie cinétique turbulente q” 


résulte de la contribution de tous les termes de I'équa- 
tion précédente, ella est en quelque sorte un constat; 
mieux vaudrait peut être considérer que le rôla des «ter- 
mes de transportr est de créer un courant d'échange 
allant des régions productrices d'énergie vers celles 
les plus défavorisées. C'est alors une relation entre 
adiffusionh et production turbulente qui serait recherchée. 
De toute maniêre toute tentative consistant à lier de façon 
plus ou moins complexe ce terms à un quelconque gra- 
diente du champ turbulent redonnera à l'êquation son 
caractêra parabolique et ce point reste à noter, La consé- 
quence d'un tel choix sur le type même de lI'équation 
étant acquise on peut dire en bref que si du point de 
vue théorique l'association d'un terme provenant du 


caractêre spécifiquament non linéaire de l['équation 
A 

Uu q avec un autre pu; dont te caractêre mixte 
CU Xy ás 


est évident" et déjá criticable, du point de vue expéri- 
mental les choses ne s'arrangent guêre mieux; en effet, 
dans une conduite Laufer*** montra que dans la région 


ou x.' est inférieur à 40 les termes uq et up sont de 
signe contraire, ce dernier étant en valeur absolue plus 


grand, ils s'annulent pour x.” = 40 puis les corréla- 


tions usp tendent à s'effacer comme l['on se rapproche 
de l'axe de la conduite**, 


Ed 


er 


Si en fin de compte, et à condition de ne pas 
séloigner par trop des cas classiques, les contrôles 
axpérimentaux demeurent tolérants, quand bien même 
une formulation hasardeuse est à la base du procédé 
cela tient gu fai; que dans le calcul des couches 
limites, jets... le rôle de ces termes n'est pas toujours 
évident. || n'en est plus de même lorsqu'ils deviennent 
la seul moteur du phénomêne étudié. Le contaminant 
«information indéfiniment fractionnable» est véhiculé 
autour de l'écoulement moyen par convection turbulente 

at 


lopérateur u, | chaque point pollué devenant à son 


UM 

tour point PA Le mécanisme est cumulatif et à 

aucun stado du processus des fractions de contami- 

nants même  faibles ne peuvent être négligées. Notons 

en passant que la portée des informations est liée à la 

puissance «trois des fluctuations de vitesse, Dans une 

couche limite ou les gradients transversaux sont de loin 
d 

les plus importants on retiendra de 3x 
ha 


U q les termes 


d o 


u, tu” + ui + us) qui comporte entre autres la 


DX 


puissance cube de us. Ainsi, alors que de part et d autre 
de la valeur moyenne, la continuité du débit, les puis- 
sances impaires d'ordre directement supérieur contrôlent 
une information qui sera par exemple dissymétrique 
dans le cas de la figure. Les petites balles fluides ani- 
mées de vitesse transversales intenses introduiront 
loins dans la direction des x; croissant une caractéris- 
tique initialement localisées en x;'. Bien entindu, le même 
mécanisme jouera dans I'autre sens mais avec une effica- 
cité moindre, Il est visible que ceci suppose un accrois- 
sement des échelles comme indiqué sur la figure. Si les 
quelques réflexions précédentes introduisent quelques 
aspects du mécanisme de diffusion elles restent plus 
critiques que constructives; sans le résoudre, la voie 
tracés par Taylor apportera d'intéressants compléments 
d'information, 


do (x) 


a Ni 


On sait que le champ de pression est liá à l'ensemble du champ de vitesse par une équation du type Poisson 


Conscients de cette diffcultá et gradant en mémoire le rôla répartiteur de la pression, rôle confirmé par la nullitá des 


corrélations pression vitesse en un point dans le cas oú la turbulento est isotrope, on a pu penser que ces corrólations se faisalent 


1 
'écho de l'anisotropia du champ turbulente exprimée par Sp ses q e prise en charge de façon convenable. 
3 


me 


H s'agit en particuliegr de letude de Laufer sur les óácoulements en canal plan et en conduite circulaira. 
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Ainsi que nous l'avons remarqué, le fait de lier la 
diffusion à un gradient donne aux équations leur carac- 
têro parabolique mais d'autres hypothêses ont été avan- 
cées entre autres par Bradshaw [1]. L'équation devient 
du type hyperbolique et il est possible de définir en 
principe une zone à I'intérieur de laquelle est localisée 
information. Il est intéressant de noter que dans cer- 
tains modêles statistiques on a pu, en introduisant une 
mémoire (donc en quelque sorte un mécanisme à relais), 
construire des équations de ce type”. 


11.2 LA THEORIE STATISTIQUE DE TAYLOR 


Nous ne rappelerons de la théorie développée 
par Taylor [34] que les aspects les plus directement 
liés à notre étude renvoyant pour de détail soit à des 
ouvrages à caractére didactiques (Tennekes-Lumley) 
soit à des articles originaux. 


Dans son étude sur la diffusion par mouvement 
continu Taylor considêre un champ turbulent station- 
naire*** dépourvu de vitesse porteuse et s'intéresse à 
travers les coordonnées de position X, d'une particule 


issue du point défini par le vecteur a 


— "t ; 
X; la, t) = a; + v; la. t)) dr 
a õ 


à des propriétés moyennes que l'on est tenté de définir 
en premiêre approximation par (X; = 8; |. On a évidem- 
ment à tout instant: 


1 d == fa 
under (X; - a; = IX, cas ã; | v; la, t) 


relation qui lie la dérivée à l'instant t de la quantité 
cherchée d'une part à la position X — a, d'autre part à 


la vitesse Vi; Par le truchement de Mp a, s introduit 
lhistoire de la particule si bien que: 


SR - 


Povo) wjlto) de 
“o 


Si l'on suppose l'existence d'une corrélation dite 
Lagrangienne telle que 
vt t)=v dt vit) 
on a en posant t=t+7 


v | A; (7) = Vj (t+ 7) vit) 


at 1 d RP 
IA, —a P= vw | Ra tr) dr 
[o] 


. 


d'ol = 4 
(x; eo im éra Vi | 0 (t == T) Rj tr) d FT 


Deux 
examinées: 


formes asymptotiques sont classiquement 


— la premiére est valable pour de petits laps de 
temps t durant lesquels le mouvement est supposé três 
fortement corrélé. R.. (7) es alors voisin de l'unité et 
Ion a: 


(X,—a P=vt 


— la seconde est applicable lorsque t est au con- 
traire três grand, On admet alors qu'il existe une échelle 
Lagrangienne fixée dês que t est supérieur à une certaine 
valeur: ainsi 


at 


Dans cet approche Taylor ne considérant que la 
dépendance entre vitesses fluctuantes d'une même par- 
ticule prise à deux instants différents, abandone le 
point de vue d'Euler pour adopter celui de Lagrange. 
Or, c'est dans l'expression de l'accélération convective 
que se situe la difficultá centrale, la présentation de 
Taylor apparait donc comme un moyen de la tourner****, 
Cependant, la mesure d'une corrélation lagrangienne 
s'avére três difficile et la possibilité d'ávaluer les cor- 
rálations lagrangiennes R à partir des classiques corréla- 
tions eulériennes R n'est guére plus simple. 


Une rapide incursion dans la domaine spectral nous 
permettra d'entrevoir le rôle des differentes structures 
dans le mécanisme de diffusion. Corrélations et spectres 
lagrangiens sont en effet liés par 


“b 4] 
E, In)=4 | o R; tr) cos 2wnrdr 
En fonction du spectre lagrangien, la résultat de 
Taylor s'écrit done: 


ia = sin” q nt 
E = =! z nn 


Ceci équivaut à dire que la dispersion à partir 
d'une source ponctuelle continue est influencée par la 


* Une chaine de Markov est dite d'ordre «ry si dans une suite d'épreuves le résultat de chaque épreuve dépend des résul 


tats des répreuves qui la précéde directement. 


**  Signalons |ouvrage récent de Monin [24]. 


+++ ne fonction aléatoira est dite stationnaire si comparant cette fonction à l'instant t avec cette même fonction à lI'instant 
t + + cette comparaison est indépendante de t. L'homogénéité est évidemment une propriété spatiale du champ. 


++** Cos corrélations lagrangiennes cherchent à définir la cohérence du mouvement d'une même particule entre deux instants. 
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totalité du spectre ce dernier étant toutefois pondéré 
par um filtre passe bas dont la frégquence de coupure 
évolue avec le temps. Ce point est essentiel car il 

nous révêle le rôle limité des frequences élevées sur 
la diffusion, or c'est précisément au niveau de ces 
structures que la reproduction en soufflerie du spectre 
atmosphérique est déficiente, compte tenu des échelles 
de réduction adoptées la viscosité rend impossible la 
simulation de structures atmosphériques dont la dimen- 
sion serait inférisure au mêtre. La remarque précédente 
montre fort heureusement que cette lacune est de peu 
d'importancce vis-à-vis du phénomêne de diffusion, 


Rermarques 


| — Dans la cadre de analogias moléculalre la dif- 
fusion avait été mise sous la forme 


aC ae õ ac 


Kia 


On peut montrer que +, est reliá à la variance de 
position des particules (X; me )º: 

Tt 

14 a 


A a; 


à titre d'information, appliquons sans précaution à ce 
cas l'analyse de Taylor qui suppose, rappelons-le, [homo- 
géênéité et la stationnarité; on trouve: 


— pour un laps de temps court 


(X—aV=v; Tr = E 


— et pour un intervalla de temps long: 
(X; E a; )* — 2 V tc; —. e = CSt 
ce qui apparenterait ce dernier cas à un processus de 
diffusion moléculaire, 


|| Le théorême central limite, dont les exten- 
sions sont du reste nombreuses, fournit une information 
sur la distribution en position des particules pour des 
temps suffisamment longs. On sait en effet que quelque 
soi: la fonction de répartition de variables indépendantes 
(| Soumises à des conditions três générales, la somme 
$=tt+ito.. +14 est assymptotiquement normale 
(n-— cv). Si Ion applique ce raisonnement au mouve- 
ment d'une particule dont la position aléatoire est défi- 
nie par: 


—— "t 
Xi la, t) =aj + | 
o“ O 


= 
v; (a, t) dt 


il faut donc tout d'abord supposer que la fonction v est 


stationnaire dotée d'échelles intégrales dliférentes de 
Zero. 


Chaque parcours résultant du découpaga peut être 
marqué à travers vi par le mécanisme turbulent, il n'en 
demeure pas moins que la densité de probabilité de la 
coordonnée de position X est pour un temps três long 
(condition assymptotique exigeant un grand nombre 
de cheminements élémentaires chacun contrôlé par une 
réalisationv, ) gaussienne. La condition assymptotique 
exige un grand nombra de réalisations contrólês chacune 


par une vitesse v. 


échelle lagrangienne (, demeu- 


t 

k-] 

rant finie de maniére à ce que l'indépendance des épreu- 
ves Ssoit respectêe. 


.3 Relations entre caractéristiques Legrangiennes 
et Eulériennes 


Les relations de base sont: 
V (a, t) = U! x (a, t) K 


m m 


: a: e” 
X ta, t)=a+ | Va ty dr 


La vitesse V de la particula issue de a est à I'instant 
t rotéa V (a, t) et sa position est désignéa par X (a, t) 
de tello sorte que V = ax (a, t)/ót. Les variables de 
Lagrange au temps t dépendent des valeurs du champ 


Eulérien U (X, t) (composantes u,) en tout point de la 
trajectoire aléatoire. Cette trajectoire est définis comme 
la solution d'un systême d'équations intégralas 


a “1 SE 
Kila, t=a, + | u 4X la, t)p dt 
Et 


ou u; est une fonction connue de la variable X. La liaison 
est en fait três complexe. Si [on admet qu'il est possibla 
d'utiliser les développements en série de Taylor on a: 


Sc a é vê | 
O: a, t) 
es ja = -— + 4 ) E 3 
Xmúl=-a+ à ——|[t-b)" 
n=10 pru t=t, 


oO = 


Via t=U la, t)+ 


|| est possible de transformer les coefficiants de 
ces séries en utilisant les relations: 


Nous admettons sans discussion que la dispersion de points mouvants est reprérable par marquage (colorant ou particules 


en suspension de densité adéquate) le mouvement de ces traceurs est identifiable au mouvement d'un point du fluide qui occupe- 


rai la position de la particule si elle ne s'*y trouvait pas. 
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q" t=t pn 
o" Yy (a, t) 
ot" t=t 
o d n Ea 
“ —>|X=a 
= 4 = dal UM dA, tt 
| E q 0X; t=t 


Lumle;y a pu tirer de ce genre de développement 
certaines conclusions, c'est ainsi qu'il montre que si le 


champ U (X, t) est statistiquement homogêne il en est 


ce même du champ vectoriel aléatoire X (a, t) mais la 
méthode tourne court dês qu'il s'agit de définir les 
caractéristiques statistiques Lagrangiennes du champ à 
partir des Eulériennes sauf pour des laps de temps três 
couris. 


En fait la plupari des auteurs introduisent qeulpues 
hypothêses simples de similitude, On peut imaginer qu'en 
identifiant un certain élément et en le suivant, la perte 
de cohérence (corrélation lagrangienne) est plus faible 
que lorsque cette même perte est observée sur deux 
éléments différents. Grosso-modo on pourrait dire que la 
corrélation eulérienne combine la décohérence du mou- 
vement de chaque particule avec une décohérence 
spaciale évaluée sur une distance qui dépend précisé- 
ment du cheminement aléatoire suivi lequel est entre 
autres lié au niveau de turbulence. Dans des cas exce- 
ptionnels tels ceux fournis par la turbulence de grille 
(ou u!/u = 2%) 

L 
pica a 
v wu 


on trouve 


On imagine donc que 
Ri tr) > Rj tr) 


|| sen suivra que léchella de temps lagrangienne 
(6 doit être plus grande que |Eulérienne T 


.- sé o 
SA + “ O 


[ =T 


On pourra également comparer les microéchelles: 


1,2 12 
Es | 2 |? 
A; E | Le Er () ) ; Ai == |- "ar (o) | 


et lã encore on s'attend à ce que A > A inégalité du 
reste confirmée par Corrsin (1963) qui précise A/ A — 
0,8 


ll ny a en fait aucune raison théorique qui permette 
d'espérer que les fonctions de corrélation Lagrangiennes 
et Euláriennes aient la même forme. Cependant les mesu- 
res récentes de Snyder et Lumley [29] laissent entre- 
voir une bonne coincidence, dans ces conditions on 
comprend qu'il soit tentant de jeter un pont entre les 
deux aspects du problême au moyen d'un simple chan- 
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gemeni d'échelle des temps comme le font Hay et 
Pasquill [14] 


RiT)=R (7) avec =Br 


Ces auteurs adoptent la valeur | — 4 
Les choses ne vont pourtant pas simplement car 
si ce résultat est confirmé par bons nombres d'auteurs 


dans le cas d'une atmosphêre calme la valeur 2 est 
êgalement avancée dans les cas instables. 


Faute de mieux, cette relation peut être appliquée 
au résultat sur la dispersion 


as a d+ 
E DE F 
(X; a; = vo t | à 


L'impossibilité d'établir une relation entra échealla 
lagrangisnne et eulérienne ne doit pas faire oubliar 
certains résultats valables pour un champ turbulent 
homogêne, le fluido étant toujours supposé incompres- 
sible. On démontre en effet que les fonctions caracté- 
ristiques et par conségquent aussi les densités de proba- 
bilité des champs de vitesse lagrangiens et eulériens 
sont identiques. Ceci implique que 

q.=m 

Toute tentative pour atteindre vi; directement devi- 

ent de ce fait inutile. 


Cependant dans un canal d'axe O x,, en admettant 
homogêne le champ lagrangien dans une section, on a: 


RE 1 o e ne 
V = ; | | uy” + (U, E U, 1º] dx: dx 
1 q r 


car alors la variance de la vitesse lagranglenne reçoit 
une contribution à la fois de la variance de la vitesse 
eulérienne et du carrá de la différence entra la vitesse 
moyenne eulérienne et la vitesse débitante U,. 


1.4 Diffuston dans un écoulement à gradient de 
vitesse constam 


Le modêle se complique alors du fait qu'il existe 
une vitesse porteuse telle que: 


VU, =A x E U=U,=O0 


Du point de vue cinématique ce type d'écoulement 
est un cas particulier de celui étudié par Craya [7] 
qui suppose À, = ô U, [dx m: La turbulence est homo- 
gêne dans le plan x; x: perpendiculaire à U,, si la station. 
narité est envisageable du point de vue énergétique 
puisque production et dissipation existent, du moins 
ella n'est pas acquise, il faut noter en effet que cet 
écoulement qui comporte une échelle de temps ne com- 
porte par contre aucuna échella de longueur autres que 
cella que l'on peut former en faisant intervenir la vis- 
cositê cinématique. On rique alors de se trouver dans 
une situation analogue à cellg rencontrée en écoulement 
derrigre grille ou toutes les échelles évoluent lentement 
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avec la distance. Nous n'en supposerons pas moins que 
la champ eulérien de vitesse est homogêne dans toutes 
les directions de sorte que le champ lagrangien est 
stationnaira. Les résultats suivants sont à prendre en 
considération ne serait-ce que du point de vue qualitatif. 


La position du particula est défini à [instant t: 


- longitudinalament par: 


a 
X, (-, |) =a + | [rx (a, vm) +v (a, a de 


NO 


— at transversalement par: 
E st 


3 * Ú 


Compte tenu des approximations faites le champ 


da vitesse v est stationnaire. Appliquant le théorême 
central limite on peut conclure directement que la distri- 
bution en probabilité des coordonnées de position trans- 
versala X; est gaussienne. Longitudinalement la situation 
est plus complexe. Considérons deux partícules issuas 
ou même instant de deux points infiniment voisins. 
Si l'une d'entre elles est écartés de la ligne da courant 
moyenne par l'action de la composant transversals vs 
tandis que l'autre suit cette ligne de courant la gradient 
de vitesse entre en jeu et de toute façon sépare les deux 
éléments par un mécanisme autrement plus efficace que 
celui résultant de action de la composante v.,. 


Poursuivant le calcul on trouve: 


RAza+ rat 


la position moyenne se déplace avec une vitesse corres- 
pondant à la vitessa moyenne relativa à la position 
initiala, 


La position relative d'une particule par rapport à 
cette position moyenne est: 


d cms 
—— (Xi — X)=A (Xi — a) +yw, 
dt 


Or, la variance de X; — as croit linéairement avec 
le temps d'aprês un résultat précédement acquis et appli- 
cable dans le plan 2 3. Aux grands temps on aura donc: 


A (Xs as) A v, 
et une nouvelle différentiation donnera: 


dº - dX. 
cms (Me Mi) À — “Ny 
dt dt "a 


Ceci montre que la quatité XxX -— "X est donnêe par 
una doubla intégration portant sur une fonction supposée 
stationaire v,, sa distribution assymptotique est donc 
gaussienne sa variance étant d'aprês Corrsin 


PE ESSE 2 
(X: = x) — CE A Ve r Css 


La dispersion dans la direction x, crolt donc beau- 
coup plus vite que dans les deux autres directions x: 


et xs pour lesquelles on a encore: 


Oy = ar = 2 vtbs 


Il est possibla de prédire la forme moyenne d'une 
tache de polluant dispersés dans un écoulement à gra- 
dient constant. La figure ci-dessous en donne une idée. 


H.5 Dispersion dans un canal 


Le canal de hauteur 2 h a son axe dirigé suivant 
Ox, de sorte que U. = U, = 0. L'écoulement moven est 
défini dans un plan x; xs. La position d'un point mobile 
est donnée par: 


= ij Pt png e 
X la =a+ | [ve 601 
a õ 
Lu EX (8, 4),4] de 


e A at 
X. fa, t] = a + | vo (a, 7) dr 
"o O 


= + —+ 
HM [a, t]=a + | va la, 7) dr 
“O 
u, désignant bien la fluctuation de la vitesse eulérianne 
à Iendroit ou se trove le point mobile. 


On peut admettre que v. est stationnaire, dans ces 
conditions X: — a: est l'intégrale d'une fonction station- 
naire (qui n'est pas elle-même obtenue par dérivation 
d'une autra fonction stationnaire), dans ces conditions 
les propriétés statistique de X. — a. sont gaussiennes. 


Les contraintes imposées par la paroi rendent la 
situation beaucoup plus complexe lorsqu'il s'agit d'étu- 
dier les propriétés statistiques de X, — as. Un raison- 
nement portant à la fois sur le caractêre stationnaire 
de X. — a; imposé par la paroi et présumé stationnaire 
dae vw qui en est la dérivée met en défaut [application 
du théorême central limite. 


Portons notre attention sur la dispersion dans le 
sens de |'éêcoulement 


X=a+U,.t 
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U, désignant toujours la vitesse débitante. Par sous- 
traction il vient: 
a EE 
W— XE | 


E [X: (a, )] — Uy + 


E e. | 
iu [X la, 7), 17] | dr 


sans entrer dans le détail disons que la stationnarité de 
X. — as fait que la point aléatoire prend les valeurs U, 
et u, suivant un processus stationnaire*. L'intégrale por- 
tani sur une fonction aléatoire stationnaire (non obte- 
nua par dérivation) les propriétés statistiques de X, =X, 
sont gaussiennes. La variance de X, — X. est donnée 
par: 


De RT = 297 Ela 


Or de 1 ME sé 1d [REO E 

Vi = 2 h | E h Wi E EU, EE U,J' dx 
Da plus 

U, — d us 

— =tfiy) et . Mu (mn) 
dk f 
Vel 
tj — x:!h 


ce qui donne: 
E | 4 “4 
vo = Uh; | É dy + | g dy | — U', 
do Q E [1] 


Il resterait à évaleur C,. Nous avons par ailleurs 
admis que (; est liée à I'échelle spatiale eulérienne L et 


à V u', Nous devons alors supposer une correspondance 
entre L et h. Dans ces conditions puisque 


TO 


L 
os U 
u 


V 


et l'on peut écrire 
(X — X)P— hUt 


Cette relation ne serait de toute façon valable que 
pour 


t>> h/U, 


Cetta suite d'exemples montra comment une com- 
plication du modéle conduit à une perte d'informations 
que seule |'expérience permet de recouvrer, Dans le cas 
de la diffusion en canal nous ne possédons que des 
formes assymptotiques, les constantes restant à ajuster, 
|| nen demeure pas moins que de l'examen de ces 
configurations se dégage les points suivants: 


— Importance dans les modéles de diffusion des 
caractéristiques lagrangiennes du champ turbu- 
lent dont les plus essentiellas s'avêrent si diffici- 
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les à reliegr à un ensemble de mesures conduites, 
à quelques exceptions prês, dans la perspective 
d'Euler, 


— Róôóle fondamental du spectre mettant en évidence 
[efficacité des basses et moyennes fréquences. 
Ceotto notion évoluant elle-méême au cours du 
temps du fait de I'importance du groupement nt. 


— Influence du grandient de vitesse qui renforce 
considérablement la dispersion longitudinala par 
Vintermédiaire des cheminements transversaux, 


— Action de la vitesse de frottement dans les écou- 
lements entre parois, 


Fauto de mieux, l'ensemble des résultats acquis 
devra donc nous servir de guide pour aborder un pro- 
blême infiniment plus complexe celui de la diffusion dans 
IV'atmosphêre. Deux nouvelles difficultés surgissent en 
afitat: 


— la premiêre résulte de l'inhomogénéité longitu- 
dinala du champ étudié qui met en échec l'appli- 
cation de certains théorêmes d'analyse établis 
pour des types bien particulior de fonctions 
aléatoires. 


— la seconde provieni des effets thermiques 
associés bien que toute notre analyse ne porte 
que sur des interactions faibles. 


|| ne faudra donc pas s étonner des approches 
semi-empiriques entre autres cell de Sutton qui postu- 
lant d'emblêe la forms de la corrélation lagrangiesnne 


écrit: 
P n 
(1) = =— 
R , p + ur 


relation dont le plus grand avantage est de demeurer 
cohérente quand on opére des passages à la limite. 
|| déduit um coefficiont de diffusion transversale qu'il 
applique à la diffusion de quantité de mouvement ce 
qui lia la répartition de la vitesse donnés sous forme de 
loi da puissance à l'iexposant n. Postulant una distribu- 
tion transversale de concentration gaussienne, Sutton 
donne une relation permettant pour une émission de 
débit Q opérée au niveau du sol de connaitre la concen- 
tration moyenne 


C (x Xs, XX) = 


20 | n—2 
== a exp. 4— x; 
me eU, x Ar | 


Ka Ma 
Eu Es Ê 


Les coefficients principaux de diffusion sont des 


fonctions de U, u, vet n, En dépit de toutes les faibles- 
ses de cet approche la relation précédente s'est déja 
avérée fort préciguse pour les applications. 


La présentation que nous donnerons ultérigurement 
des résultats obtenus s'appuiera sur la notion de simili- 


*On voit bien le rôle important joué par les parois qui confinent |écoulement. 
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tude lagrangienne. Sur le plan matériel nous quitterons 
en rêgle générale les notations indicielles pour revenir 
aux appellations classiques x =x (direction du vent), 
x = j (direction transversalo), x» =z (direction verti- 
cala) la concentration moyenne étant toujours désignée 


par € 


Batchelor admet que dans la zone turbulente interne 
d'une couche limite, lã ou la contrainte tangentielle 
n'est pas trop éloignée de T ps les propriétés statistiques 
du mouvement d'une particulas émise au niveau de la 
paroi peuvent être décrites en ne faisant intervenir que 
Upettet bien entendu z; si la paroi est rugueuse ce qui 
est géênéralement la cas dans las problêmes qui nous 
occupent. Pour un écoulement stratifiê Gifford [12] puis 
Cermak [2] introduisent comme paramêtre supplémen- 
taire la longueur de Monin Obukhov. Pour un grand 
nombre de réalisations la position X Y Z d'une particule 
au bout du temps «th est telle que: 


dX U, z Zn | 
merge ei E ms 1 = 18 pi ) 
k L L 


b et c sont des constantes universelles introduites par 
Batchelor qui propose b = 0,1 e e — 1, Les fonctions f et 
y devraient être universelles et reliées au profil de vites- 
se moyenne. 


De nombreux auteurs ont par la suite proposé des 
extensions ou des modifications de cette formulation. 
En particulier, à la suite d'un calcul d'Ellison (1959) 
recoupant des mesures de Cermak (1963), il est admis 
que b — k; par ailleurs, Monin (1953) a établi une 
importante relation à partir de considération d'énergie: 


1 | 1/4 
ts = | — —— 
| 


Compte tenu de ces remarques, la théoris de la 
similitude lagrangienne peut founir une évaluation quan- 
titativo dês que I'on se donne une forme de profil de 
vitesse, c'estã-dire dês que f est connu, Il est alors 
possibla de calculer la trajectoire du centre de masse 
d'un nuage de particules marquées, et la loi de décrois- 
sance de la concentration maximale au sol. Ainsi, Pano- 
fsky et Prasad (1965) ont remplacé f par une évaluation 
numérique, tandis que par exemple Gifford (1962) utilise 
la loi de Kazanskii et Monin: 
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dans laquelle 
no: cal | = 0.03 
L 


Enfin, plus récemment, Pasquill (1966) considere 
que les résultats de cette théorie doivent mieux s'appli- 
quer aux cas réels si, au ligu de Z, la cote moyenne 
du panache dans une section est exprimée par Z telle 
que: 

C (x, O; £p] 


Es =-0,1 
Cmax (x;) 


1 


£, est ici supposé mieux rendre compte que Z de la 
forms de la distribution verticale de concentration, 


CHAPITRE Ill 


DIFFUSION DANS UNE COUCHE LIMITE STRATIFIEE 
RESULTATS EXPERIMENTAUX 


Hl.1 Conditions d'essais 


Nous allons rappeler três briêvement dans quelles 
conditions se sont déroulés les essais de facon à définir 
les états atmosphériques auxquels ils correspondent et 
à mieux faire comprendre I'intérêt de dispositifs prévus 
au départ pour élargir la gamme des modêles simulés 
(injection par paroi poreuse loin en amont). 


L'installation est représentée sur la figure 1, nous 
ne nous attarderons ni sur sa description, ni sur celle 
des méthodes de mesure les détails les plus importanta 
étant donnés par Solal [31]. Nous ne rappelerons que 
las points suivants: 


— deux vitesses d'écoulement ont été expérimentées 
UV, =4,3m/s., U/*= 6,6 m/s. 


(la soufflerie autorise 2m/s<U, ?P < 20m/s. 
“la longueur de Monln Obukhov L variant de 


—- % à — O,4m 
(la soufflerie autorise: 
U,, 
é = ess do UM) Ei 
k — Hp 
T p Cp 


La vitesse de frottement est calculêe par la for- 
mule de Ludwieg et Tillmann, la paroi inférieure 
équipée en plaque fluxmétrique permet de défi- 
nir Hp. La température T, est prise égale à la 
température moyenne dans la couche, 


- L'épaisseur de la couche limite évolua cans la 
région da mesure autour d'une valuer dc évaluée 
à 15 cm. Si l'on estime à 30% la fraction de la 
couche ou la tension reste à peu prês constante 
et êgale à 7 " la partie intéressante est de I'ordre 
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de 5 cm. Des mesures ont été effectuées jusqu à 
des distances x, , pouvant atteindre 2,5 cm en 
aval du point d'injection; le développement du 
panache dans le sein même de la couche se pré- 
sente alors comme indiqué sur la figure ci-contre, 
|| ne semble cependant pas que la zone turbu- 
lente externe ait été appelée à jouer un rôle 
déterminant. Cette remarque n'appelle donc rien 
d'autre qu'une certaine prudence dans l'inter- 
prétation des résultats correspondant aux sec- 
tions les plus éloignées du point d'injection, elle 
fait mieux comprendre l'intérêt porté par maints 
auteurs à des procédés cherchant artificiellement 
à épaissir avec la zone turbulente pariátale. Bien 
qu'il n'ait pas été utilisé lors de ces essais, qui 
s'inscrivent en fait comme une toute premiêre 
tentative la soufflerie ici décrite en comporte 
un dont nous décrirons par la suita quelques 
caractéristiques essentielles. 


iu 
4% e Fe 


| 
[o - e % 
| .% % AA Polluze 


Po url dource 


— Le gaz traceur utilisé est de I'éthane, les concan- 
trations sont détectées au moyen d'un analyseur 
opérant par ionisation de flamme, L'injection se 
fait aussl discrôtement que possible au niveau de 
la paroi. 


Conditions de similitude réalisées 


Nous nous proposons d 'analyser biéêvement à quel- 
les conditions atmosphériques réellas correspondent les 
mesures entreprises. Rappelons tout d'abord qu'une 
bonne récapitulation des donnés concernant l['atmosphêre 
a été donnée par Slade (1968). 


On peut se faire aisément une idée de I'échelle gêo- 
métriques adoptée. Elle se fonde essentiellement sur les 
dimensions respectives de la zone à tension constante. 
Remarquons qu'elle occupe dans latmosphêre environ 
50 mº. Dans notre cas oú |'éápaisseur totale de la couche 
limite est 8c = 15 cm la même région s'étend vertica- 
ment sur 0,35c = 5 cm. Le facteur d'échelle géométri- 
que adopté est donc 1/1000 environ. 


*La valsur adoptée pouvant aussi atteindre 150 m. 
we de 


Ce facteur multiplicatif doit être pris en compte pour 
toutes les grandeurs caractéristiques De ce fait il devient 
difficile de simuler en souffleris des écoulements forte- 
tement instables. L'écoulement le plus instable que nous 
avons réalisé jusqu'àã présent correspond à L-—- — 0,4 m, 
il reproduit donc une instabilitê atmosphériqua corres- 
pondant à L = — 400 m et par conséquent três faible 
puisqu'une atmosphêre caractérisé par une valuer L = 
— 100 m n'est encore qualifis que de «légêrement 
instable», l'appelation instable à três instable correspon- 
dantãá ILI < 10 m, Notons qu'il est possible d'atteindre 
avec les moyens actuels L = — 0,1 m ce qui correspon- 
drait dans l'atmosphêre à . |LI| = 100m, 


En ce qui concerne [échelle de la turbulence on 
peut affirmer que, eu égard aux dimensions de la couche 
limite, les strutures turbulentes dont l'échella caracté- 
ristique est supérieure à 1 m ne véhiculent probablement 
que três peu d'énergie, Ceci signifia que la dimension 
des plus grandes structures atmosphériques qu'il est 
possible de simuler sont de l'ordre de 1000 m. Si l'on 
compte avec une vitesse de vent correspondante égale 
à 5 m/sec. le temps caractéristiqueo que Ion peut alors 
formé est de l'odre de 200 secondes. La question da la 
similitude spectrale entre atmosphêre et souffleria se 
pose alors avec acuitéá sachant par ailleus le rôle joué 
par cette répartition énergétique sur le phénomêna de 
diffusion. 


Examinons la courbe spectrale relevée dans l'atmos- 
phére par Van der Hoven (1957) Fig. 2. Par chance une 
large bande de fréquences comprise entre quelques 
heures et quelques minutes ne contient que três peu 
d'énergie. Cette remarque capitala va permettre en quel- 
que sorte de découper le phánomêne en deux parties qui 
selon toute vraissemblance ne s'interactionnent que fai- 
blement I'une l'autre**. C'est uniquement la partie haute 
frequence de ce spectre mouvementé que l'on s'efforcera 
de réaliser en souffleria, La fréquence de coupure que 
Ion peut accepté correspondent à 10 cycles/heure 
(temps caractéristique égal à 6 mn) n'est que le double 
de cells que la souffleria décrite permet de reproduire 
en l'absence de tout artifice. Au paragraphe suivant nous 
verrons de façon plus quantitative comment se présente 
en fait cette simulation spectrale. 


Apreés avoir abordé les problêmes relevant du domai- 
ne des basses iréquences examinons les difficultés que 
Ion pourrait rencontrer à l'autre extrémité du spectre 
dans le domaine des hautes fréguences, Si l'on suppose 
dans la couche turbulente pariétale (zone à 7 = 7 cons- 
tant) 


— d'une part un équilibre entre production et dis- 
sipation, 


— d'autre part l'isotropie des plus fines structures 
on a: 


La continuité des spectres relevés en turbulenue ne peimet pas un tel découpage. Si une telle situation se prósentait on 


pourrait aucueillir avec infiniment moins de suspicion |'hypothêre d'Heisenberg qui tente fictivement de l'opérer, les tourbillons cor- 


respondant à des nombres d'onde supérieur à | k | se comportant vis-á-vis de ceux dont le même nombre est inférieur comme une 


viscosité fictive. 
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o OU U ' 15 vu” 
Wo =UW— = — et y, = — 
' dz kz A” 
d'oú U'; À 
Wo = —— = cste 
zru 


On peut s'attendre à ce que 1 u*/ U, dépende três 
peu du nombre de Reynolds, le mécanisme turbulent étant 
du point de vue production contrôlé par U, par contre 
si I'on désigne par U* une vitesse de référence le 
rapport U/U* en dépend. Cette connexion reste toutefois 
assez faible, négligeons-la pour l'instant. Dans ces con- 
ditions l'expression précédente donne: 


UT A A 2 2 
] z Re= cste avec Ri, = - 
A) Z [y 
Ainsi 
j E 
À N FR EM 


| | am a 


z soufflerie z aimosphêre 


Du fait que la similitude de Reynolds n'est pas 
respectée ceci entraine 


A A 


e E 


z soufflerie z atmosphêre 


Dans le cas qui nous occupe, si Ion admet que la 
vitesse de référence est reproduite dans la maquette 
(par exemple U = 5 m/sec, ce qui correspond à un vent 
de l'ordre de 20 km/h) le fluide utilisé étant air dans 
les deux cas le rapport des nombres de Reynolds devient 
celui des échelles géométriques soit 1 000. Le rapport 
a/z sera environ 39 fois plus petit en soufflerie que dans 
l'atmosphêre; lors de comparaisons ultérieures on obser- 
vera cette différentiation des spectres pour les fréqueces 
élevées mais Il'on sait que ceci est sans conséquence 
grave pour les problêmes de diffusion. Ces considéra- 
tions sur la similitude appeleront une dernigre remarque 
relative aux temps d'intégration pratiqués en soufflerie 
et dans l'atmosphêre. Nous avons en effet précisé que 
les temps caractéristiques des strutures turbulentes 


était en soufflerie de éme de seconde ce qui corres- 


pondait dans l'atmosphêre à 3 mn environ. Or, la plupart 
des mesures qui y sont effecluées portent sur des temps 
d'intégration de l“ordre de I'heure. Nombre d'auteurs se 
sont inquiétés des problêmes posés par une telle trans- 
position. On peut en dégager les formules approchées 
suivantes que nous utiliserons dans notre cas: 

C max [3 min] 2 [3 min.) GE 


ese RO =p = 
C max [1 heure] o [1 heure] 


S pura A A 
2/1 heure] - 


Ces différents examens montre en fin de compte 
que si la simulation en soufflerie demeure possible elle 
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appelle bien des précautions et Ion voit que les essais 
qui vont être présentés se situent à la limite des réalisa- 
tions fiables. Il est certain que |'épaississement contrôlé 
de la couche aura une influence extrômement bénéfique, 
la construction et l'implantation de souffleries de plus 
grande dimension posant quelques difficultés pour des 
équipes ce recherche à caractére purement universitaire, 
la longueur de plancher de celle utilisée avoisine 10 m, 
sa section étant supérieure à 1 mí. 


Amélioration des procédés de simulation 


|| n'est pas dans notre intention de les décrire tous 
tant ils sont nombreux et par bien des points décevants, 
Lorsque l'on désire obtenir une répartition de vitesse 
ayant un profil domné il vient immédiatement à I'idée 
d'utiliser un grillage à porosité variable dont la déffinition 
ne présente aucune difficulté. Malheureusement si le 
champ movyen est conforme au but à atteindre le champ 
turbulent ne l'est pas, L'espoir de voir ce dernier rapi- 
dement se restructurer en présence du champ de vitesse 
de façon convenable est vite decu ce qui semble en 
partie condamner le procédé, du moins sous sa forme 
élémentaire. J. Mathieu propose un soufflage à travers 
une paroi poreuse située loin en amont des sections 
de mesure et le résultat semble satisfaisant ainsi que 
nous allons le voir à partir des résultats de Schon et 
Mery [30] au moins tant que le soufflage reste modéré; 
nous ne chercherons en fait d-ns ce qui suit qu'ã com- 
parer des résultats, obtenus avec et sans soufflage, 
l'épaississement recherché ne dépassant pas 100%. 


Du point de vue des vitesse repésentées figure 3 
dans un systême de variables convenables, la coinci- 
dence est três satisfaisante. Les intensités turbulentes 
longitudinales et transversales n'accusent que des écarts 
faibles du moins dans la région turbulente pariétale 
qui nous intéresse et à condition de se trouver à plus 
de 6 m en aval de la surface d'injection (Fig, 4). La 
situation parait tout aussi favorable si I'on compare la 
répartition des contraintes de Reynolds avec ou sans 
injection; (Fig. 5) une réalisation convenable de cette 
caractéristique est bien évidemment cruciale puisqu'elle 
conditionne avec le profil de vitesses movennes les 
valeurs du coefficient de diffusion e. 


Les spectres de ul, ví et w' ont été mesurés.Nous 
donnons figure 6 les deux plus importants dans une 
représentation en n F [n] qui renforce le rôle des struc- 
tures moyennes dont l'importance est considérable dans 
la mécanisme de la diffusion, La coincidence avec les 
spectres relevés dans l'atmosphêre est tout à fait remar- 
quable sauf pour les fréquences élevées ce qui est 
conforme à notre attente, La situation n'est pas plus 


critique pour des valeurs de nz/U = 0,005 ce qui, compte 
tenu des échelles de transposition adoptées, correspond 
dans l'atmosphêre à quelques cycles par heure. Dans la 
représentation choisie la zone inertielle du spectre corres- 
pond à une pente — */, elle est plus étendue comme 
il fallait s'y attendre pour le spectre en u” que pour 
celui en «w, 
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| faudrait toutefois souligner qu'en dépit de ces 
répartitions énergétiques satisfaisantes les intensités 
turbulentes convenables ne sont pas respectées la paroi 
étant lisse, 


Bref, parmi les moyens d'épaissir artificiellement 
la couche limite tout en conservant les caractêres essen- 
tiels propres au problême de diffusion celui du soufilage 
amonit est à considérer avec grand intérêt si Ion se 
contente de doubler au plus l'épaisseur de la couche. 


[11.2 Examen des résultats expérimentaux 


Les résultats seront présentés sous une forme 
classique qui présuppose une structure voisine de I'équi- 
libre. 


Dans une section donnée repérée par son abscisse 
x, lorigine coincidant avec le point d'injection, las 
échelles caractéristiques sont au nombra de trois: 


— |a concentration maximum dans le section Cmax 


— deux échelles de longueur? et respectivement 
définies par 


€ e 2245/2,52) 


=:01 
Cmax [x] 


Les coefficients numériques ont été choisis de 
telles façon que si la concentration est gaussienne 9 
représante |'écart type. 


La distribution de concentration est assez bien 
représentés par 


RE ER 
= exp. (= 043 (5) “080 (5) 
E 


[o] 


Cmax 


comme le montre la figure 7. Une loi normale correspon- 
drait évidemmant à la relation 


IC Ra: ZA E 
= exp -os0(5)—os0( 5) 


Cmax E, z 


L'évolution des échelles longitudinales Cmax, 2 et 
cá reste à définir. Du fait que le phénoméne semble 
en gros refléter un état d équilibre on la cherchera 
d'emblée sous forme de loi de puissance adoptant une 
triple représentation logarithmique. La concentration 


Cmax est introduite sous la forme adimensionnelle 
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Cmax U, 
8] 


ouú O est le débit injecté. Les figures 8, 9 et 10 mettent 
en évidence l'influence sur la concentration de |'éloi- 
gnement et de l'instabilité, L'évolution est assez bien 
représentés par des relations du type 


Crmax -no. 


los trois exposants p, h et v dépendant de L. Sur? q et 
“a on note même une légêre influence de la vitesse 
laquelle semble traduire lI'influance du nombre de 
Reynolds ce qui n'est guêre surprenant du fait que la 
paroi sur laquella se développe la couche est lisse et 


Vinjection faite au niveau même du plancher. 
Dans le cas d'une distribution gaussienne un con- 
tróle du débit injecté conduit à la relation suivante 


entre les exposants 


p=h+(m+ 1) v-— ma 


Ici on a pris a 
5 cin. 


U (2) E 


3 
TE: | et cin. — ae 


De ce point de vua três global les résultats présen- 
tés paraissent corrects les faibles écarts semblant davan- 
tago imputables à des erreurs accidentelles qu'ã des 
erreurs systêématiques 


L=-— 0» L=—-04 m 
p calculé 1,48 1,53 
p expérimental | 1,42 1,52 


|| est possibla de cet ensemble de résultat de tirer 
une évaluation de es. On conduira le calcul en se plaçant 
dans le cas d'une linéaire équivalente à la source ponc- 
tuella donnêe 


+ 00 
CM 2) dy 
source ponctuella 
— Pe 


C sourcelinéairo (x, Z) = | 


Entre la paroi et la cote z la débit de polluant 
convectéê est: 


“Z 
p (kw £) = | C (x 2) U bu 2) de 
“O 


cequi fixe la débit de matiêre «,,s échappant vertica- 
lement 
d 
To = p (xi Z) 


dx 


On peut alors écrire: 


Ad .Z a 
E | E (% 2) U fã, 2) de 
CMT v6) 


E ix, 2) = a 


e (xs, Z| 


dz 


La Figure 11 fait état des résultats trouvés et tas 
compare aux formes proposées. Le rôle destabilisant 
du chauffage apparait comme capital. 


Déásirant comparer plus en détail entre autres ce 
dernier résultat avec la méthode de calcul mise au 
point au laboratoire par Jeandel et Mathieu nous indi- 
querons seulement que des confrontations ont été faites 
par Solal [31] entre ses résultats et la méthode de 
calcul préconisée par Sutton. L'accord est satisfaisant. 


Nous avons donné l'essentiel de la théorie de la 
similituda lagrangienne, On peut en retenir que cette 
formulation permet de calculer de façon três simple la 
concentration maximale à partir de la seuls connaissance 
trajectoire moyenne d'un ensemble de particules et la 
du profil de vitesse. Les résultats sont présentés sous 
la forme adimensionnelle préconisés par cette théorie 
figure 12 et figure 13. Pour la soufflerie utilisée la hau- 
teur de rugosité a été évaluée fictivement à z, = 10“m. 


Les courbes en traits pleins sont des estimations 
faites par Gifford (1962) et Panofsky-Prasad (1965) 
à partir de diverses évaluations du profil de vitesse dans 
la couche limite atmosphérique. On peut noter la bonne 
concordance entre les valeurs pour l'atmosphêre et les 
valeurs mesurées en souffleria. 


La méthoda de calcul que nous avons développés 
a pour point de départ une idés émisa par Kovasznay 
et Neo [15] selon laquelle la coefficient de diffusion 
turbulente peut du point de vue de son évolution être 
assimilé à toute quantité turbulento qui est convectés, 
diffusée, produite et dissipée. Si l'on tient compte 
du fait que la cosfficiant « s'apparente à celui des corré- 
lations uo en un point, il est vraissemblable qu'il puisse 
vérifiar un tel schéma, Un des avantages de cette formu- 
lation est sans doute d'éviter une évaluation toujours 
contestable du terme diffusion car dans le modéle initial 
proposé « contrôlait, à la convection prês, l'ensemble du 
processus qui assure la transmission de l'information à 
travers la couche., Il est intéressant de constater que 
dans cette perspective on retrouvait la caractêre non 
linéaire du problême lã même oú nous l'avions localisé,. 
En bref |'éguation proposés prenait la forme 


tre — de o dve 9U, 

=— 4 y mi A Que =v) —— 

Us dx, Us 0X OX * dg + ' Tt 
— o e me e 
Convection par le Diffusion Production 


mouvement moyen 


Bye (vo -v) 


a 1 
L d 
O gi 


Dissisation 


OU Po = à + p 
Le terme production de », était non sans analogie 
dUi 
U; 207 ASSOCIÓ à es, et à un gradient 
E] 


de vitesse; un coefficiente multiplicatif A devant pe-- 
mettre un ajustement aux données expérimentales. 


avec la terme u; 


Quant au teme dissipatif s'inspirant des lois de 
décroissance de la turbulence isotrope et remarquant 
qu'en l'absence de tout autre mécanisme (production, 
conduction, diffusion) on aurait aussi 


= — diseipation vc 
ou vc est introduit sous la forme 
Ps ==9] Pá 


ou appa:ait une nouvelle constante B également à 
ajuster. Quant à la longueur L, en principe liée aux 
échelles dissipatives, elle sera supposée variar comme 
x; prês de la paroi et comparable à &c dans la région 
turbulente extérieure. Bien entendu cette approche sous 
entend une proportionnalité entre échelles macroscopi- 
ques et microscopiques. 


La méthode révéla bien vite certaines faiblesses 
en particulier lorsque d'importants gradients de pression 
entraient en jeu, les profils de vitesse ne vérifiaient pas 
les données expé-imentales et la croissance de la couche 
limite était trop lente. La possibilité de prendre en 
charge l'ensemble de la diffusion au sens large (enten- 
dons par lã lensemble du mécanisme assurant la trans- 
mission des informations à travers la couche, la convec- 
tion par la mouvement moyen exceptée) parut alors 
suspecte. || existe dans une couche limite de vastes 
mouvements qui sembleraient devoir être traités à part. 
Le cheminement des informations par la biais de ces 
structures à grande échelle relêverait alors davantage 
d'un mécanisme convectif que d'un processus de con- 
tamination plus ou moins contrôlé par des paramêtres 
locaux. L'adjonction d'un terme 


Poxa Nº / rr, max. dpi 
Ae | = a U, o e 
| 80 UP do Ox 


fut donc envisagé. Ce terme fait nécessairement inter- 
venir le gradient de »., son efficacité restant bien 
entendu liés à ['inhomogénéité transversale du champ 


* Méthodes da résolution des áquatisas aux dérivêss partielles de la couche limite. C. R. Acad. Sc. t. 274, mai 1972: 
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turbulent dve/dx,. L'intensité de ce terme à caractére 
plus ou moins convectif est ajusté à la vitesse de réfé- 
rence U, * tandis qu'une double pondération est opérée. 
La nécessité de voir s'effacer avec r, tout processus 
turbulent appelle son introduction comme facteur multi- 
plicatif par exemple sous la forme adimensionnelle 
“;max./U 9 8c. La seconde pondération fait intervenir 
la distance x;/ãc à la paroi. ÁAinsi complétée I'équation 
se met sous la forme: 


à [| dr U B (», — 4] 
ae [e O 

Le choix des constantes et de la fonction L| a été 
opéré comme suit: (Fig. 14) 


Are = 8,2.10%; A = 0,133; B-0,8; L = 


La contribution comparée des deux termes assurant 
dans une section la transmission de [information par 
voie turbulente est donnée sur la figure 15. On peut 
s étonne: de voir croitre avec z/ãc un terme qui d'une 


cetaine façon relêve de la fluctuation «w' de la vitesse; 
or Ion sait que cette quantité décroit lorsque I'on s'éloi- 
gne de la paroi. Cependant si l'on tient compte des 
fonctions spectrales relatives à chaque point et ne s'inté- 
resse qu'àã la fraction d'énergie correspondant aux struc- 
tures turbulentes dont les nombres d'onde restent com- 
pris entre O et 27/5c on s'aperçoit qu'il s'agit bien lã 
d'une fraction croissant de z/5c ce qui sans justifier 
Vévolution du terme la rend vraissemblable (Fig. 16). 


En l'absence de stratification, la répartition de ec 
trouvês par cette méthode était en trés bon accord avec 
les mesures de Solal (Fig. 17). 


Pour tenir compte de l'effet destabilisant du chauf- 
fage pariétal nous avons dans un premier temps admis 
que dans l'équation précédente seul le terme production 
était affecté. Un coefficient multiplicatif «ks fut intro- 
duir en relation soit avec le nombre de Richardson soit 
avec la longueur de Monin. Dans l'un des cas taités 


(L= — 0,4m) nous avons pris k = 1,5 dans |I'équation 
ainsi modifiée: 
; Z 4M2 (vo. max. d 
UE +wW— + ae ( - e 
dx dz 8 US 3c iz 


FE dro! dU (pe— ») 
o Ba o) + M (pe — p) 2 — B O Pe 


Compte tenu de ces aménagements un três bon 
accord existe entre la répartition calculás et mesurée 
de «e». Ici du reste dans cette premiêre approche nous 


* Dans le cas présent wj sera nul, 
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n'avons pas tenu compte des écarts existant entre coef- 
ficient de diffusion de quantité de mouvement et coef- 
ficient de diffusion de masse lesquels peuvent dans 
les cas choisis différer de 10% en rágima neutre et 
de 20% en regime faiblement instable. 


La méthode de résolution du systême 


ou dW 
do — Ea) 
dx dz 
ou dW d du 1 aP 
U — = ve | — — — 
dx dz dz | “ dz | p dx 
et 
i dp dp. oz 1 o 
U ada fa er 
dz dz | e dz 
d Oro ' du tre Es: v) 
ss do o ps a 
x R 2 A k Veia v) E B T am 


s inspire beaucoup de celle préásentée par Spalding et 
ses élêves [32]. L'introduction d'un maillage t;ans- 
versale en relation avec la fonction de courant assure 
sa constante adaptation à la dimension de la couchea. 
La variable 


| ] 
= 4 


O) — 


er, 
E a] 


se p.ésente comme une fonction de courant adimension- 
nelle ouy et + figurent les valeurs limites externes et 
internes* Avec comme nouvelles variables indépendantes 
x et w le systême précédent prend la forme: 


4)? do 


DU amp OU O [ PU», du 1 daP 
e “vai dx 


ou mp désigne le taux de transfe;it de masse à travers 
la frontiêre externe de la couche. 


Ox ç E* v | LE = d | du 
| pé U ve dre p ou 
q e A kt,—r)— E 
do |lve— yr | - VE-4I du 
DME 
LU 


Il est alors intéressant de remarquer que les deux 
équations précédentes prennent la forme commune 
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dep dep d dp o 
E a o — = + da 
dx te fa ping + | s dus ) 


Us du 


ou & représente soit U soit »,, a est un teime qui varie 
suivant que I'équation se rapporte à U ou »,: il en est 
de même de ab considéé comme terme «sources. 
|| représente dans I'équation de quantité de mouvement 
action des forces de pression dont |'évolution suivant 
Ox est fixée par I'écoulement sain; dans l'équation de 
comportement de », il traduit |'áquilibre turbulent pro- 
duction dissipation. De toute évidence les premiers ter- 
mes prennent en charge la convection longitudinale et 
transversale ainsi que les transports tu-bulents à grande 
échelle, le troisiême tient compte du processus de diffu- 


sion transversale. 


Cette similitude de forme donnée aux équations 
du mouvement est três intéressante sur le plan du trai- 
temente numérique. Bien entendu, une étude particuliere 
de la rágion pariéátale doit àtre fait le modele ne devenant 
convenable que lorsque z/&c est supérieur à 0,05. Les 
fonctions algébiiques qui sont employées par Spalding 
at Patankar (1967) et qui sont fondées sur une modifi- 
cation de Ihypothêse de Van Driest ont été reprises. 
Les détails du calcul sont donnés pa: Jeandel* et seront 
publiés ultérigurement avec ensemble des résultats 
acquis par Jeandel et Mathieu. 


CONCLUSION 


Le parti le plus intéressant que l'on puisse tire: de 
cet ensemble de réflexions et de résultats est, nous 
semble-t-il, de les comparer avec ceux obtenus dans 
Vatmosphêre. Outre le problême des temps d'intégration 
choisis qui diffêrent énormément d'un cas à l'autre, 
même si l'on fait jouer la notion d'échelle (3 mn p.r. à 1 
heure) une correspondance entre les paramêtres emplo- 
yés doit encore intervenir. 


La tu.bulence atmosphérique n'est pas dans les 
mesures pratiques définies avec un fil chaud. L'instru- 
ment utilisê est assez souvent une girouette à faible 
temps de réponse qui enregistre simultanément les fluc- 
tuations des deux angles 4, etá. repérant la direction 
horizontale et verticale du vent. Les grandeurs mesurées 
sont alors les écarts types cor:espondant 46 et de: De 
plus6p joue souvent le rôle d'indicateur grossier de la 
stabilitê aérothermique de |'écoulement. Une estimation 
de la valeur desgcorrespondant, en moyenne, à I'écou- 
lement reproduit en soufflerie a été publiée à l'aide d'un 
tableau par Cramer, De Santo, Dumbauld (1964). Ce 


tableau relie la valeur de 4g, à celle des exposants «h» 
et «v» des lois d'évolution en fonction de x, des échelles 
transversale 4 et verticale 4 


Fe fi k [ 
| cs 
4 0,56 0,86 
5 0,64 0,88 
6 0,71 0,91 


L'accord entre les résultats expérimentaux trouvás 
et les valeurs sg et 64 déduites de mesures effectuées 
dans l'atmosphêre est três encourageant. (Fig. 18). 


A lissue de ce premier bilan puels sont les avan- 
tagens comparés des mesu.es en soufilerie et des mesu- 
res «in situs. Nous n'avons manqué de signaler les diffi- 
cultés d'une transposition provenant entre autre du fait 
que les échelles de temps réalisées en soufflerie demeu- 
raient três faibles. La couche à tension constante s avê- 
rait exiguê pour une soufflerie de la dimension de celle 
utilisée. Enfin et surtout les instabilités réalisges pou- 
vaient être considérées comme peu significatives bien 
que nous ayons déjá pu constater leur influence consi- 
dérable sur le coefficient de diffusion turbulente +. Il est 
certain que toutes ces objections sont fondées, cependant 
il n'est pas difficile à la lumigre des tentatives précé- 
dentes de percevoir les voies susceptibles de conduire 
à de notables améliorations. Sans modifier la géométrie 
du tunnel il est facile, par soufflage, de doubler I'épais- 
seu; de la couche; en jouant à la fois su; cet artífice 
et sur la longueur du plancher on peut espérer multiplier 
par quatre cette dimension, Toute réduction de la vitesse 
est intéressante si l'on che:che à accroitre l'instabilité 
de lécoulement. On est bien sôr limité dans cette voie 
pa: les difficultés métrologiques diverses qui sont d'im- 
portance surtout si l'on veut connaitre la structure fine 
de Iécoulement par exemple à travers les corrélations 
vitesses températu:ie. Une surchauffe de la paroi, bien 
qu'envisageable, reste au moins dans notre cas limitée, 
ella demeure de toute façon bien moins efficace qu'une 
action sur la vitesse qui si l'on néglige l'effet du nombre 
de Reynolds intervient à la puissance trois dans la 
définition de L. En conjuguant les deux moyens il doit 
être possible de réduire d'un facteur 4 la valeur de L. 
En conclusion, les commodités relatives d'une expári- 
mentation qui atteind déjã avec précision les flux et 
les concentrations et qui nous permettra três prochai- 
nement de connaitre certains aspects de la structure 
turbulente de l'écoulement ne doivent pas être perdu 
de vue. Quand bien même ce type de recherche ne 
saurait se substituer aux mesures effectuées dans l'atmos- 
phére il semble bien qu'il soit le plus apte à déceler 
des tendances et à contrôler des modéles. 


* Méáthode de calcul de la couche limits turbulente. Thése Lyon 1972. 
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FIG. T — Coubes adimensionnelles de distribution de concentra- 
ton vorticale et Iransversale 
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FIG. 11 — Coefficiant da diffusion verticale de masse 
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Analysis of rain rate distributions in Portugal 


and its application to the microwave 


propagation above 10 ghz. 


RESUMO 


Este trabalho representa mais uma tentativa para se 
obter o conhecimento dos regimes de precipitação instan- 
tânea em Portugal, com vista ao cálculo de ligações em 
freguências superiores a TO GHz. 


A análise estatistica apresentada foi conduzida 
com base em registos de precipitação cedidos pelo 
serviço Meteorológico Nacional. Tais registos foram 
obtidos com pluvidmetros do tipo de sifão com uma 
precisão que permite leituras em intervalos de tempo 
de 5 minutos. 


1 — INTRODUCTION 


The growing applications of the microwave band 
above 10 GHz, has determined the need of propagation 
measurements under rain conditions. 

Also Po tugal develops a research program under 
the auspices of Junta Nacional de Investigação Cien- 
tífica e Tecnológica which intends to study the various 
effects of the rainfall on the propagation at such fre- 
quencies. 

These p:opagation measurements will be carried 
out at Ovar and Lisbon. 

While installing both stations it is important to 
know: 


fi) The attenuation to be expected. 
(iil The possibility of extrapolating the data obtai- 
ned at these stations to other places. 


This work is an attempt to answer to both ques- 
tions. For this purpose we tried to find out some correla- 
tions of the type established by other authors, some of 
which have already been tested (1) a few months ago. 

We shall divide this report in three main parts: 


(il Analysis of the registrations available at the 
Portuguese Meteorological Office, in order to 
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SUMMARY 


This work represents an attempt to estimate rain 
rate distributions at different points in Portugal. These 
distributions are of great importance in order to calculate 
the attenuation to be expected on paths working above 
IO GHz. The analysis of the rainrate frequencies was 
based upon the registratinos available at the Portuguese 
Meteorological Office, obtained by means of syphon 
type rain gauges, which allow to measure the total am- 
ount of precipitation during 5 minutes intervals. 


obtain the 5 minutes cumulativa rain rate dis- 
tributions. 


fii) Correlation of the results mentioned here before, 
with various meteorological parameters. 

(iii) Calculation of the attenuation distributions using 
the results presented by other authors. 


2 — ANALYSIS OF THE RAIN RATE REGISTRATIONS 


This analysis was carried out with a 5 minutes 
integration time, as mentioned above. 

In order to reduce the registrations to be analised, 
we assumed that only the 5 minutes rain rates above 
24 mm/hour were significant. This means a total amount 
of 2 mm during the 5 minutes interval. 

These procedures where applied to the 13 mete- 
orological stations of the Portuguese network where 
it was possible to obtain registrations with an accuracy 
compatiblo with the 5 minutes measurements. 

Fig. 1 to 4 show the mean rain rate cumulativa 
distributions for a period of ten years from 1961 to 1970. 

Fig. 5 and 6 represent the correspondent distri- 
butions for two stations whose registrations only cove- 
red a period of 5 years. 
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3 — CORRELATION OF THE RAIN RATE DISTRIBU- 
TIONS WITH VARIOUS METEOROLOGICAL PARA- 
METERS 


As mentioned above we choosed 13 stations of the 
Portuguese network where it was possible to obtain 5 
minutes rain rates. However the Portuguese net work has 
about 200 stations, where various meteorological para- 
meters are obtained. 

Some of them can be correlated with rain rate fre- 
quencies as some authors have already done. 

Among all the parameters available at each station 
we choosed: 


— Annual amount of precipitation. 

— Number of days with measurable precipitation. 

— Number of days with a total amount of precipi- 
tation higher than 10 mm. 


—— Number of days with thunderstorms. 


First of all we tried to establish the regression lines 
of average daily precipitation for days with measurable 
rain versus the number of 5 minutes inte.vals that rain 
rate exceeds the specified intensities. 


Fig. 7 shows this result for the 13 stations. 5 minu- 
tes intervals are represented as a percentage of the year. 

We conclude from fig. 7 that these parameters are 
not correlated. As a matter of fact, the selected para- 
meters do not individualize the stations, as we find the 
same average daily precipitation for some sites with 
different regimen of precipitation. Then, we represented 
the number of days with thunderstorms versus the num- 
ber of days with a total amount of precipitation higher 
than 10 mm. 

Fig. 8 shows the correlation between the two men- 
tioned parameters established for all the stations of the 
Portuguese network where they were available, 

We can conclude that there is a higher dispersion, 
mainly for the stations of: 


Porto (Pilar and S. Gens) 
Coimbra 

Braga 

Vidago 

Caramulo 

Moncorvo 


Viana do Alentejo 


Finnaly we tried to correlate the cumulative distri- 
butions of 5 minutes rain rate directly with the number 
of days with thunderstorms and also with the number 
of days with a total amount of precipitation above 
10 mm. 


it was not possible to establish any correlation. 
However there are some important conclusions to be 
pointed out: 


(i) The upper limit of the annual highest 5 minu- 
tes rain rate lies between 115 and 125 mm/h 
except for the Caramulo Station where it is 
150 mm/h. 


(ii) The lower limit of the annual highest 5 minutes 
rain rate lies between 40 and 65 mm/h. Once 
more the Caramulo Station is an exception and 
the lower limit is 90 mm/h. 


(iii) The cumulative distributions of rain rate are 
very different from station to station and for the 
same station in different veas. 


4 — CALCULATION OF THE CUMULATIVE DISTRIBU- 
TIONS OF THE ATTENUATION 


As it is shown in (1) it seems that the distributions 
given by Boithias et al (2) allow us to calculate the 
cumiulative distributions of equivalent rain rate along paths 
with different lengths. So, using the distributions pro- 
posed by Boithias et al combined with the Ryde's law 
and with the average distributions of rain rate measured 
at the mentioned 13 stations, we calculated the average 
cumulative dist:ibutions of the attenuation at 11 GHz for 
two paths located in the area of those stations. 


Fig. 9 to 21 show the distributions that we obtaincd. 
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Sistemas de segunda ordem 


resposta no tempo e na frequência 


RESUMO 


Algumas experiências simples permitem aos alunos 
dos cursos de engenharia electrotécnica visualizar as 
relações entre a análise matemática e o comportamento 
fisico de sistemas lineares. 

Demonstra-se, a partir de um circuito RLC série, as 
relações entre a resposta no tempo e na frequência. 
Sugerem-se vários tipos de experiências laboratoriais 
de possivel realização pelos alunos. 


1 — INTRODUÇÃO 


À equação diferencial de segunda ordem com coefi- 
ciantes constantes aparece frequentemente em electro- 
técnia. O seu estudo torna-se muito fácil se nos debru- 
çarmos sobre o circuito RLC série 


Ro. L 


FIG. 1 


À equação diferencial deste circuito é 


do db 


MO=IRAL + 


dt (1) 


a as suas soluções foram ao longo do curso apresenta- 
das em vários cadeias (*) com especial incidência 
para respostas a uma tensão escalão e uma tensão 
alternada sinusoidal. 


2 — RESPOSTA DO CIRCUITO RLC SÉRIE A UM 
ESCALÃO 


Consideremos o circuite da figura 1 e apliquemos 
a entrada o escalão de tensão V 1 (t) em que 
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RICARDO CHARTERS D'AZEVEDO 


SYNOPSIS 


Simple experiments engineering students to realize 
the relationship between mathematical analysis and the 
physical behavior of linear systems. 

The teorical and experimental evidence is used to 
present the relationship between time-domain and fre- 
quency - domain response. 

Suitably simple laboratory experiments are provide 
for the students. 


144) 
1 (t)=1 t>0 (2) 


t<O 


Com para t= O não circula corrente, podemos 
escrever 


V CV 
A A RA AGA O LOGO RG 1 (3) 
onde 
| lp) é a transformada de Laplace de | (t) 
L [i(t)]=1(p) com t variável real e 


p variável complexa ou simbólica. 


A corrente tem uma lei de variação no tempo 
diferente conforme os valores da resitência em relação 
à capacidade e à inductância. As raizes do trinómio 
do denominador podem ser reais ou complexas depen- 
dendo do binómio 


A=R'C'—-4 LC (4) 


24—A>71 lió. R>2 VL/C — regime sobre- 
periódico) 


Às raizes são reais e da forma: 


Í / L 
=-a=-— À 144/ 1-4—— 
l EL  Y RC (5) 


* Este trabalho tem sido realizado por alunos do 4.º ano de Engenharia Electrotécnica. 


a 


R / L 
= Bs: — 1 — 1 oo: — 
; dia 4 V RC (6) 
e a corrente terá a forma 
V 
p= === 
L (p+a) (op + b) (7) 
com 
R L 
Ea Eee quão pe / 1-4 — 
É: L RC (8) 
A transformada inversa é 
—bt — at 
V 
o e == ei 4 — 6 
pm 
/ 4L 
R / 1 so — YE 
V'- qe (o 


A tensão aos terminais do condensador é dada por 


1 1 “tê 
V — | E A s i(t) dt 
cte) o (0) a | (10) 
ou seja 
—bt — at 
ce — be 
Val)=V (1 — mmpeeviiarmnia 
ft a-b ) (11) 


A tensão parte de zero com uma tangente horizontal 
para atingir um valor limite igual a V, sendo a carga 
tanto mais rápida quanto menor fôr a resistência (sendo 


lógicamente sempre superior a 2 VL/C). 


22—4=0 (ié. R=2VL/C- regime aperió- 
dico limite) 


As raizes são iguais e dadas por 


R 
'=D0'=-gãa-— —m— 
p p 2L 
e a corrente é dada por 
V 
| — 
(p) Na 
E es 
LL (12) 
À transformada inversa dá-nos 
— Rt/2 L 
V 
(tg) = — te 
L (13) 
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e para a tensão aos terminais do condensador vem 


= R/Zl 


R 
4 t—— 8 


ve 8 2L (14) 


Vet) =V 


cor espondendo à carga mais rápida sem ultrapassar 
o valor da amplitude do escalão. 


23—sc<Olié. R< 2VL/C, regime periódico) 


Escrevemos a coorrente | (p) sob a forma 


V 1 
da RN2 71 Rº 
per — + |—— — ) 
og À LE ar. (15) 
8 como 
1 Rº 
LC” qu 
fazemos 
1 Rº R 
—ems — o à =——— 
LS 4 LU 2 
e obtemos 
V — at 
i (t) = sen ot 
Ta E (16) 


Se o termo expononclal e igual à unidade, ou seja 


mi 
R<< 2 V c não há amortecimento e a corrente sinu- 


soidal circula indefinidamente com uma pulsação 


Quando o amortecimento não é desprezável, a cor- 
rente mantem um andamento sinusoidal, mas cada má- 
ximo tem uma amplitude menor: que a precedente. 

Define-se então uma «pseudo-oscilação» 


=V TETO (17) 


8 uma constante de amortecimento 


R 
QU sao 
ak (18) 


Para os instantes 


qu pp tn =3, ds (cai) 
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a corrente passa por zero mas os seus valo:es máximos, 
obtidos pela anulação da derivada de it), não estão q = 
situados exactamente em 7/2u; 3m/2w... se bem que 
afastados uns dos outros dum número inteiro de meios 
períodos (T/2 = m/w) — fig.2. 


R (21) 


Neste caso a carga do condensador é oscilante e O 
valor V de amplitude é atingido ao fim de um certo 
número de passagens positivas e negativas em rela- 


cão a V. 
— qt 
E sen ut = v) 
NV) T+ =— RT 
| s=ng (22) 
[1] 
Com tag = — 


Os instantes para os quais a tensão é igual ao 
valor V, são dados por 


t— p=nT 


FIG. 2 
ou seja 
A razão entre a corrente no instante t, e o valor 
da corrente no instante t. é dado por q 1 td 
=: aims arc tg — 
ti pr] dd 
— tl 
Is o sen ut; a, o que corresponde a valores ligeiramente superiores a 
e: Sonia =: 
R pa art. m/2u, 37 f Bús.. 
e sen ut: (19) 
| Quanto aos valores máximos são obtidos pela anu- 
pois que lação da derivada de v (t). Notemos que o valor máximo 
wt.= ut + 27 ni A 
Vmde =>V(1+e ““máx]) 
donde máx 
log-— =0eT=- =) =V(1+exp-—— pese 
tu 2 11 (20) á ' (23) 


Define-se assim o decremento logaritmico à co pode atingir 2V se a constante de amortecimento a 
obtem-se o factor de qualidade Q fôr nula. 


3 — RESPOSTA EM FREQUÊNCIA 
Se considerarmos a função saída/entrada obtemos 


v5 1 


— (p= ———— 
o CO 8,4 B6+41 (24) 


ou de uma maneira mais geral 


r5 tdo 
ve adia p'+HZep+ mw (25) 
Considerando o número y = — Vo , que pode 


tomar todos os valores de zero a bm dependendo 
dos valores da resistência em relação à indutância e 
capacidado teremos, 


FIG 3 31—y>1 (R>2YL/C) 
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Podemos decompor o trinómio do denominador 
num produto de dois termos 


rs 1 


re e: p 

mos a a 

a | b ) (26) 
em que a e b tem os valores dados em (5) e (6). 


Como a > b obtemos as assintotas representadas 
na fig. 4 


FIG. 4 


Para R >> 2VL/C obtemos 


R 1 
Eu — ES —— 


L RC 


1 R 
«a . Então para vw = 1/RC 
RE “LL ESA 


obtem-se uma atenuação de cerca de 3 dB. 


muito afastados pois 


32 4 =4 (A= 2/6) 


FIG. 5 
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O denominador de (25) é neste caso um quadrado 
perfeito. 


le 1) (27) 


Obtem-se então duas assintotas que se cruzam no 
ponto (O dB; w = wo) 


A frequência de corte a — 3 dB é definida para 
ILCo+1)=2 
ou LC w = 0,643 


0,643 
ou fo | 
2x LC (28) 
Sendo por definição o tempo de crescimento, o 
tempo que ve lt) leva a evoluir entre 1/10 e 9/10 do 
valor final obtemos para este caso 


tc = 0,34 LC (29) 


e obtemos 


Neste caso as raizes são imaginárias puras 


p= ju 
, 8 vem 
LS Pa 
sa ET 
O = À cão” TE vw”) + 2 Jjaw (30) 
-ib (wu) 
=A lol » (31) 


Designando-se por À (w) a característica de resposta 
em amplitude e o seu argumento — b (w) a característica 
de resposta em fase. 


Considerando 
a lub=-— log A lu) 
e exprimindo em dB vem 


rã 
A ah 


— à (0) =10 log —-—————————— 
na (og — 0)" + 4 au 


e a constante de resposta de fase virá dada por 


Z2euw 
b (wo) =arc tg — - 
Vo — w 
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Quando a frequência tende para zero a (w) tende 
também para zero dB, e quando a frequência tende para 
infinito a (w) varia com — 40 log LCw; obtemos assim 
duas assintotas 


—— uma recta horizontal de ordenada O dB 


— uma recta a — 40 dB/década (ou-12 dB/oitava) 


Reparemos, no entanto, que a curva real se afasta 
tanto mais das assintotas quanto ; se aproxima de zero 


Aus) 


Ru 


9,54 (<07307 


“dal. 


FIG. 6 


Podemos determinar alguns pontos notáveis para 
uma caracterização mais perfeita da resposta em ampli- 
tude. 


Para w = uol 1—-2% a derivada é nula e a (w) 
passa por um máximo para y < 0,707 


a (w) dB é nulo para 


o = O 


us = us V2 A = 9º) 
Soluções só possiveis se y < 0,707. 


À curva de alo4p corta a assintota 12dB/oitava em 


ia 
ú = 


V2 (1—242) 


também para y < 0,707. 
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1 
Podemos reparar QUE «way ws CÊ : as frequências são 


simétricas em relação a « (em coordenadas log) e a 
frequência limite superior — 3 dB é definida por 


E À be == 


! 


= am (1— 2%) V + 


Paray <. 1 a resposta inicial apresenta sempre uma 
sobrelevação mesmo quando *; está compreendido entrre 
0.707 0 1. 


4 — TRABALHO A REALIZAR 
4.1 — Dimensionamento prévio. 
4.1.1 — Para os valores dos parâmetros represen- 


tados na figura, calcule os valores de R que levam aos 
seguintes valores de y = 0,3; 0,4; 0,5; 0,6; 0,7; 1; 2; 5. 


FIG. 7 


4.1.2 — Determine os valores das raizes no plano 
complexo p e a sua localização para cada um dos valores 
de 7 determinados atrás. 


4.1.3 — Calcule a máxima tensão que pode surgir 
aos terminais do condensador em função da tensão Ve 
de entrada. 


4.1.4 — Calcula os tempos de crescimento previsi- 
veis para os valores de 7 >1 


4.2 — Determinações experimentais 


4.2.1 — Utilizando um gerador de varrimento, um 
registador XY e um fasimetro trace as características 
de resposta em aplitude e em fase para os diversos 
valores de R calculados em 4.1.1. 


4.2.2 — Com um gerador de impulsos, um oscilos- 
cópio e a respectiva camara fotográfica regista as res- 
postas em frequência e os respectivos tc para os mesmos 
valores de R. 
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4.3 — Conclusões 


Tire as conclusões que lhe sugerem os resultados 
obtidos. Anote e lustifique as discrepâncias entre os 
valores obtidos em 4.1 e os de 4.2, 


NOTA FINAL 


A partir do estudo feito é possivel obter um sem 
número de aplicações sujeitas a prática laboratorial. 


Se se dispuser, por exemplo, de um registador de 
quadro móvel, este trabalho é realizado com um número 
muito menor de equipamento e com uma ligação muito 
mais perfeita com a teoria dos aparelhos de quadro 
móvel [1]. 


Para que este trabalho seja completo sugere-se 
que sejam investigadas as soluções a partir de um com- 
putador analógico. 


CONCLUSÕES 


Este artigo descreve um trabalho experimenta! que 
pode levar os alunos a participar alguns aspectos da 
resposta de um sistema linear. Faz-se uma introcução 
teórica no domínio das respostas na frequência e no 
tempo e sugerem-se algumas aplicações possíveis. 


Este trabalho é normalmente realizado pelos alunos 
do curso de Engenharia Electrotécnica na Academia 
Militar quando da sua passagem pela cadeira de Medi- 
das Eléctricas. 
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a 


à 
dd 
bi F 
E) sê Ê 
= 


Aa 
» a transmissão 
a corrente. * 


| O CORRENTES DE PRECISAOÚ 


di monte a transmissão 
exe wsva | completa utilizando 
OS PRODUTOS 


O ACIONAMENTOS PADRONIZADOS 


até 4250 HP à 300 RPM 
PARA A TRANSMISSÃO DE FORÇA 


| O CORRENTES 'NON-LUB' 
o ou redutores, corrente transportadora, luva de acoplamento, 


impraticavel 


O CORRENTES SIMPLES DE 
ROLOS EM AÇO INOXIDAVEL 
para aplicação especializada 


para locais de lubrificação 
embreagem, freio, variador de ça etc. 


Infe ormaçõe s tecnicas sobre à linha completa de produto ENOLD na 


HARKER SUMNER & CIA.LDA. 
so = 38 Rua de Ceuta, 48 PORTO 
| E Apartado No. 162. Tél: 27054 (4 linhas) 
SALES DIVISION 14 Largo do Corpo Santo, 18 LISBOA-2 


MANCHESTER - ENGLAND à Tél; 324823, 39124. 
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Les ressources, la production et la demande 
Mundiales D'uranium (*) 


La demande d'uranium correspondant aux divers et aux Ressources Raisonnablement Estimées (R.R.E.) 
programmes consacrés à l'énergia nucléaire dans le (3), qui sont ellesmêmes subdivisées entre deux caté- 
monde, qui s'éléve actuellement à environ 18.000 gories de prix: (a) moins de 10 dollars par livre d U,O, 


tonnes de métal par an, atteindra vraisemblablement (fmoins de 26 dollars par kg d'uranium métallique) et 
quelque 60.000 tonnes par an en 1980 et 100.000 (b) de 10 dollars à 15 dollars par livro d'U,O, (de 26 


tonnes par an en 1985. «ll est done primordia| d'adopter dollars à 39 dollars par kg d'uranium métallique) (4) 
incessamment des mesures pour accélérer les travaux Les chifffres ainsi fournis — en date du mois de janvier 
d'exploration de l'uranim, de mariêre à pouvoir main- 1973 — s'établissent, pour la tranche de prix infárieure 
tenir des réserves assurées adéquates». à 10 dollars par livre d'U,O,, à 866.020 tonnes d'uranium 


(1.126.000 tonnes courtes d'U,;0,) pour ce qui est des 
C'est lã l'un des point principaux qui se dégage du R. R. A. et 916.000 tonnes d'uranium (1.129.000 tonnes 


nouveau rapport intitulé «Uranium: Resssources, produe- courtes d'U,0.) dans le cas des R. R. E.; pour la tranche 
tion et demande», qui a été préparé conjointement par de prix comprise entre 10 et 15 dollars par livre d'U,O,, 
des Groupes de travail spécialisés de I'Agence de FOCDE ils sont de 680.000 tonnes d'uranium (884.00 tonnes 
pour |'Energis Nucléaire (AEN) et de |'Agence Interna- d'U,0.,) en ce qui concerne les R. R. A. auxquels s'ajou- 
tionale de I'Energis Atomique (AIEA). En effet, malgré tent 632.000 tonnes d'uranium (821.000 tonnes courtes 
la três rapide expansion que Ion prévoit pour la demande d'U.0.]) correspondent aux R.R.A.. 
d'uranium, le rapport précise que «les prix actueles de 
Vuranium ne sont pas, en général, suffisants pour per- Ces chiffres font apparaitre, dans la tranche de 
mettre d'entreprendre |'exploration requise ou donner prix inférieure à 10 dollars par livre d'U,O., des augmen- 
lFaccélération nécessaire à l'expansion des movens de tations de quelqua 220.000 tonnes d'uranium pour les 
production. || conviendrait par conséquet de trouver des R.R. A. et de quelque 240.000 tonnes d'uranium pour 
méthodes garantissant que I'on pourra atteindre les nive- les R. R. E.. Dans la tranche de prix comprise entre 10 et 
aux de production requis por éviter uns instabilité du 15 dollars par livre d'U,O,, on constate une augmentation 
marchá au cours des années 80”. de quelque 100.000 tonnes d'uranium par rapport aux 
chifíres fournis dans le précédent rapport (1970) en ce 

En ce qui concerne les ressources, les chiffres qui concerne tant les R. R. À, que les R. R. E. Une répar- 
donnés correspondent, comme dans le précédent rapport, tition des chiffres par pays est donnée dans le tableau 1 
aux Ressources Raisonnablamente Assurées (R. R. A.) [2] ci-joint, 

(*) «Uranium —- Ressources, Production et Demandes OCDE, Paris 1973. 


(2) Définies en tant que euranium se trouvant dans des gisements connus dont la quantitó et la configuration sont tellas 
qu'on peut, dans les limites de prix donndes les récupórer de façon rentable grâce aux techniques d'extractlon et de traitement 
actuellement éprouvées». 


(3] Définies en tant que «uranium dont on suppose la présenca dans des prolongements inexplorés de gisements connus, 
ou dans des gisements nom découverts, situés dans des districis ou l'on sait pouvoir trouver de [uranium; on présume en outre 
qu'il est possible d'exploiter cet uranium de maniére économique, à un prix situé dans la tranche considárdies. 


(4) Par suite des changements de parité qui ont affectá le dollar américain (la valeur du dollar américain en date de 
mars 1973 a été adoptée pour l'ensemble du nouveau rapport), les valeurs absolues de ces tranches de prix sont en movenne 
d'environ 17% inférieures à celles du précédent rapport. 
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quelque 370.000 GWe (370.000 kWe) pour la zona euro- 


dont 280 GWe se trouvant dans les 


pays de la CEE); 540 GWe en Amérique du Nord (dont 


500 aux Etats-Unis) 


Ê 
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| sont 


tallée, qui 
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e 100 GWe au Japon, 


Une autre question qui dépend du rythme d'expan- 
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Quant à la demande, les estimations fournies dans 


le nouveau rapport sont basées sur les derniêres pré- 
visions disponibles (avril 1973) en ce qui concerne 


représentées sur le graphique 1 ci-joint. Ces prévisions 
montrent que, vers les années 1990, la capacité instal- 


[expansion de la puissance nucl 


L'importance de ces taux d expansion, du point 
de vue de la demands de combustible, est mise en 
évidence dans le graphique 2 ci-joint ou les courbes 
A et B représentent les limites supérieures et inférieures 
que devrait atteindre la demande, la courbe À supposant 
la maintien des tendances actuelles dans le choix des 
filigres de réacteurs des (c est-à-dire celle des réacteurs 
thermiques avec une forte prédominance de ceux à eau 
légêre) et la courbe B correspondant à une demande 
légêrement inférieure, qui résulterait de [introduction des 
surrégénérateurs rapides dont le rythme commencerait à 
saccélérer au début des années 1980, Les courbes A? 
et Bi de ce même graphique font apparaitre les taux 
supérieurs e infériours auxquels on pourrait s'attendre, 
en ce qui concerne la demande de combustible, dans 
la cas d'expansions de la capacité nucléaire installée 


qui aboutiraient respectivement à une capacité totale de 
20% supérieure et de 25% inférieuro à celle prévue 
en 1990, 


Un autre point qui se dégage du graphique 2 permet 
de conclure que les variations de la demande dues à 
d'éventuelles variations de l'expansion énergétique serai- 
ent trés largement supérieures aux modifications qui 
pourraient résulter de l'adoption de stratégies différentes 
en matiêre de réacteurs. 


Enfin, le graphique 2 comporte également les cour- 
bes À, B, A2, B1, rapprochées dans le temps de 8 ans, 
pour représenter les réserves considérées comme néces- 
saires pour garantir 8 années d'exploitation et que Vin- 
dustris miniére de l'uranium s efforce de maintenir pour 
des raisons d'ordre technique et économique. 


RELATION ENTRE LES RESERVES D'URANIUM 
ET LES BESOINS CUMULES EN URANIUM 


10º tonnes d'U 
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1500 | a ss 
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* Ces courbes sont avancées de huit ans. 
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Une autre question qui dépend du rythme d'expan- 
sion de la capacité installêe et des stratégies en matiêre 
de filibres de réacteurs, est celle de la demande relative 
à [uranium enrichi at à la capacité des installations 
d'enrichissement. A cet égard, le rapport montre que 
si, en ce qui concerne lenrichissement, la capacité 
mondiale actuelle ou prévue au cours de la prochaine 
décennie semble devoir suffire jusqu'au début des années 
89 ou même jusqu'en 1985, il faudra ensuite disposer 


d'une capacité supplémentaire qui, lorsqu'elle commen- 
cera à simposer, devra vraisemblablement subir une 
três rapide expansion (voir graphique 3 ci-joint). C'est 
lã une conclusion qui revêt une três grande importance 
si [on tient compte des «edélais» (qui sont à peu prês 
du même ordre que ceux requis pour les réacteurs 
eux-mêmes) à prévoir pour la construction et la mise en 


exploitation des nolvelles installations d enrichissement. 


POTENTIEL MONDIAL DE TRAVAIL DE SEPARATION 
RESULTANT DE LA PRODUCTION DES ETATS-UNIS ET D'EUROPE” 
(Teneur résiduelle en U 235 de 0,275 4 avec recyclage du plutonium) 


193 tonmes UTS 
50 


- 200 


1974 197% 1978 1980 


— 50 


4 -100 


— 150 


— 200 


1982 1984 1986 1988 199 


Les chifires relatifs à la capacitá annoncés par Eurodif na sont pas inclus. 


EXPANSION DE LA PUISSANCE NUCLEAIRE 
INSTALLEE DANS LE MONDE * 
(Chiffres nets) 


GWe 
ss | 
x chiftres correspondants tirés du rappot de 1970 sé 
1000 Monde” O 22100 
E 
500 | 500 
200 | 200 
100 | 100 
50 50 
20 20 
x 
E AE À Ee | E a N 
1972 1974 1976 1978 1980 1982 1984 1986 1958 1990 


Dans le présent rapport, par «Monde» on entend tous les pays, hormis I'URSS, I'Europe orientale et la Chine. 
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O consórcio Brown Boveri 
na distribuição de energia 


O disjuntor é um O nosso programa de Planificamos e cons- 
elemento base na fornecimento compre- 'truimos instalações 
distribuição de energia ende: exteriores e interiores 
eléctrica -— Disjuntores pneu- para todas as tensões 
máticos até 765 kV 
- Disjuntores de peque- 


no volume de óleo 
-— Disjuntores SF, 


Os disjuntores BBC Condicionamentos Com prazer resolve- 
respondem a todas técnicos, de serviço e mos o vosso problema 
as exigências de económicos estão na 

serviço base da escolha de um 


dos três sistemas de 
extinção do arco 
eléctrico 


ia 


E 


Disjuntores pneumá- E Disjuntores de  Distunorss SF. 
ticos pequeno volume de 
elíTo) 


Sociedade de Electricidade Brown Boveri, Lda. E pe CG 


Rua de Sá da Bandeira, 481 — 2º Dto. - PORTO BROWN BOVERI 
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J. MATHIEU 
DIFUSÃO NAS CAMADAS LIMITES ESTRATIFICADAS 
Técnica N.º 425 — | —= 4.1975, p. 237 a 262 


usando uma atmosfera padrão em referência à qual é anali- 
sado o escoamento real sob a forma de uma perturbação. 
Estabelecem-se os pertis de velocidade e temperatura médias 
na camada limite atmosférica para várias condições de 
estabilidade. 
Faz-se o estudo da difusão turbulenta pela analogia mole- 
cular e pela teoria estatística sendo analisadas as dificul- 
dades decorrrentes da mão linearidade do fenómeno. Mostra- 
-S8 como uma perturbação térmica pode acelerar o processo 
da difusão turbulenta ou pode bloqueá-lo consoante o tipo 
dae estabilidade da atmosfera. 
Referem-sa os ensaios efectuados em túnel aerodinâmico, as 
limitações e vantagens em relação às medidas «in situa. 
às dificuldades de rransposição dos resultados do túnel 
para a atmosfera são devidas em primeiro lugar à pequenez 
das escalas de tempo e comprimento no túnel e à dificul- 
dade em simular regimens com instabilidade significativa. 
Apresentam-se vias para melhorar estes aspectos. 
à comodidada e economia das medições feitas em tunel 
justificam a sua utilização e desenvolvimento mesmo que 
| sejam aplicadas apenas para o estudo de tendências ou 
efeitos globais. 


| 
Estuda-se a estrutura das baixas camadas da atmosfera 
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ANALISE NOS REGIMES DE PRECIPITAÇÃO INSTANTANEA 
| EM PORTUGAL E SUA APLICAÇÃO NA PROPAGAÇÃO 
| DE MICROONDAS ACIMA DE 10 GHz. 
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Esto trabalho representa mais uma tentativa para se obter 
o conhecimento dos regimes da precipitação instantânea 
em Portugal, com vista ao cálculo de ligações em frequên- 
cias superiores a 10 GHz. 

A análisa estatística apresentada foi conduzida com base 
em registos de precipitação cedidos pelo serviço Meteoro- 
lógico Nacional. Tais registos foram obtidos com pluviá- 
metros do tipo de sifão com uma precisão que permito 
leituras em intervalos de tempo de 5 minutos 
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SISTEMAS DE SEGUNDA ORDEM 
RESPOSTA NO TEMPO E NA FREQUÊNCIA 
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Algumas experiências simples permitem aos alunos dos 
cursos de engenharia electrotécnica visualizar as relações 
entro a análiso matemática a o comportamento físico de 
sistemas lineares. 

Demonstra-se, a partir de um circuito RLC série, as relações 
entre a resposta no tempo e na frequência. Sugerem-se 
vários tipos de experiâncias laboratoriais de possivel reali 
zação pelos alunos. 
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SECOND-ORDER SYSTEM: TIME-DOMAIN 
AND FREQUENCY-DOMAIN RESPONSE 
Tácnica N.º 425 — L — 4.1975, 0. 277 a 282 
Simple experiments enginneering students to realize the rela- 
rlonship between mathematical analysis and the physical 
behavior of linear systems. 
The teorical and expsrimenta! evidence is used to present 
tha relationship between time-domain and frequency-domain 
response, 
Suitably simple laboratory experiments are provide for the 
students. 
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Jd. MATHIEU 
DIFFUSION IN STRÁTIED BOUNDARY LAYERS 
Técnica N.º 425 — |! — 4.1975, p. 237 a 262 


The structure of the atmospheric boundary layer is studied 
by means of a «deviation flow» from the standard atmos- 
here. 

Mag velocity and temperature profiles are stablished for 
several stability conditions. 

Turbulent diffusion is studied both through a gradient trans- 
por: approach and throunh a statistical theory. Difficulties 
due to the non linearity of the nhenomena are analysed. It is 
shown how a destabilizing thermal action can sneed up 
turbulent diffusion. 

Measurements in wind tunnel and its limitations and advan- 
tages are comparad with «in situ» measurements. Difficulties 
intransposig the wind tunnel measurements to full scale are 
mainly due to the small wind tunnel time and lenath scales 
and to the difficulty of simulating fairly unstable flows. 
Several ways of improving the measurements are presented. 
Economy and convenience of wind tunnel measurements 
justify its utilization and development even if our purpose is 
to seek tendencies or global effects. 
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ANALYSIS OF RAIN RATE DISTRIBUTIONS IN PORTUGAL 
AND ITS APPLICATION TO THE MICROWAVE 
PROPAGATION ABOVE 10 GHZ. 


Fécnica N.º 425 — L-— 4.1975, p. 263 a 276 


This work representa an attempt to estimate rain rate dis- 
iributions a: different points in Portugal. These distributions 
aro of great importance in order to calculate the attenuation 
to be expected on paths working above 10 GHz. The 
analysis of the rainrate frequencies was based upon the 
registrations avallabla at the Portuguese Meteorological 
Office, obtained by means of syphon type rain gauges, 
which allow to measure the total amount of precipitation 
during 5 minutes intervals. 
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